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Abstract 

We introduce simple models of genetic regulatory networks and we proceed to the mathematical analysis of 
their dynamics. The models are discrete time dynamical systems generated by piecewise affine contracting 
mappings whose variables represent gene expression levels. When compared to other models of regulatory 
networks, these models have an additional parameter which is identified as quantifying interaction delays. 
In spite of their simplicity, their dynamics presents a rich variety of behaviours. This phenomenology is 
not limited to piecewise affine model but extends to smooth nonlinear discrete time models of regulatory 
networks. 

In a first step, our analysis concerns general properties of networks on arbitrary graphs (characterisation of 
the attractor, symbolic dynamics, Lyapunov stability, structural stability, symmetries, etc). In a second step, 
focus is made on simple circuits for which the attractor and its changes with parameters are described. In the 
negative circuit of 2 genes, a thorough study is presented which concern stable (quasi-)periodic oscillations 
governed by rotations on the unit circle - with a rotation number depending continuously and monotonically 
on threshold parameters. These regular oscillations exist in negative circuits with arbitrary number of genes 
where they are most likely to be observed in genetic systems with non-negligible delay effects. 



1 Introduction 



With genome sequencing becoming a widespread procedure, the structural information contained in the 
genome is now largely accessible. Not much can be said about the functional information contained in 
gene expression regulatory mechanisms which is still largely unravelled. Under a simplifying point of view, 
regulatory mechanisms are described in terms of networks of basic process in competition. Insights on the role 
of regulatory networks during the development of an organism and changes with environmental parameters 
can be gained from the analysis of crude dynamical models |27j . The models are usually supported by directed 
graphs where the nodes represent genes (or their products) and where the arrows represent interactions 
between genes. 

Depending on context, various formalisms have been used to model regulatory networks, see jS] for a recent 
review. Discrete variable models (boolean networks) have been employed to obtain essential features, such as 
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influence of the sign of a circuit on multi-stationarity or homeostasis [2E1 EE] or decomposition into regulatory 
modules |2E]- In addition to the analysis of behaviours for arbitrary graphs, models of specific regulatory 
mechanisms have been derived and thoroughly investigated in this formalism, e.g. flower morphogenesis in 
Arabidopsis thaliana |19| and dorso- ventral patterning in Drosophila melanogaster |23| . 

In a more traditional framework, (systems of coupled) nonlinear ordinary differential equations with in- 
teractions represented by sigmoid functions and early piecewise affine analogues have been considered 
Piecewise affine models admit analytical investigation without affecting most features of the dynamics. These 
models have been studied by using tools inspired from the theory of dynamical systems. In [7j , a method 
has been developed to determine existence and stability of periodic trajectories with prescribed qualitative 
behaviour. In [H] and in symbolic dynamics has been employed in a computational framework in order to 
obtain results on qualitative behaviours and their changes with parameters, namely bifurcations. Naturally 
coupled differential equations have not only been analysed in their own but they have also been applied to 
represent specific mechanisms, see |2Z] for a review. 

In spite of being governed by the same rules (described below) , boolean networks and coupled ordinary 
differential equations often present distinct dynamical behaviours |26| . Boolean networks reveal periodic 
orbits where differential equations only possess stationary points. 

These distinct behaviours however can be recovered in a unique and simple model by adjusting a parameter, 
say a. The model, an original discrete time dynamical system with continuous variables, obeys the same 
rules as in previous models. Gene product concentrations (expression levels) evolve according to combined 
interactions from other genes in the network. The interactions are given by step functions which express that 
a gene acts on another gene, or becomes inactive, only when its product concentration exceeds a threshold. 

Based on comparisons with (systems of coupled) delay differential equations, the parameter a appears to be 
related with a time delay. When this parameter is such that the delay is maximum, the system reduces to a 
boolean network. And in the limit of a system without delay, the dynamics is as for a differential equation. It 
is therefore of particular interest to describe the phase portraits not only depending on interaction parameters 
but depending on the delay parameter a. 

That a discrete time dynamical system provides a simple substitute to a delay differential equation is a 
well-known fact which has already been employed in mathematical biology IT5] . Note that this does not 
contradicts the statements on equivalence between discrete time dynamical systems and ordinary differential 
equations (Poincare section). Indeed, such statements match the dynamics of discrete time dynamical 
systems in an A-dimensional phase space with the dynamics of an ordinary differential equations in an 
N + 1-dimensional phase space. 

The paper is organised as follows. With the definition of the model provided (section |2J), we proceed to a 
comprehensive analysis of its dynamics with emphasis on changes with parameters. The analysis begins with 
the study of general properties of networks on arbitrary graphs and more specifically on circuits (section EJ. 
In particular, attractors are described in terms of symbolic dynamics by means of an admissibility condition 
on symbolic sequences. Lyapunov stability, structural stability and symmetries of orbits in the attractor are 
presented. A special section discusses relationships with boolean networks and differential equations. 

In a second step, we focus on simplest feedback circuits with 1 and 2 genes (section 0J. An analysis of the 
dynamics is presented which is complete in phase space and in parameter space for most of these circuits. 
The results rely on previous results on piecewise affine and contracting rotations. 

The dynamics of negative circuit with two genes requires special attention and a fully original analysis. The 
most important orbits (the most likely to be seen in numerical experiments with systematic prospection 
in phase space) are the so-called regular orbits. These are stable (quasi-)periodic oscillations governed 
by rotations on the unit circle composed of arcs associated with atoms in phase space. The associated 
characteristics (rotation number and arc lengths) depend on systems parameters in agreement with the 
parameters dependent sojourn times in atoms of such oscillations. 

For the sake of clarity, the paper is decomposed into two parts. In part A (sections |2] and 2J>, we only 
present results. Some results can be obtained easily and their proof are left to the reader. Most original 
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results however require elaborated mathematical proofs and calculations. These are postponed to part B 
(sections (HJ to[7J. 

The final section |S| contains concluding remarks and open problems. In particular, it indicates how the 
regular orbit analysis developed in the negative circuit of two genes extends naturally to circuits with 
arbitrary number of genes. 



2 The model 

Basic models of genetic regulatory mechanisms are networks of interacting genes where each gene is submitted 
both to a self-degradation and to interactions from other genes. The self-degradation has constant rate and 
the interactions are linear combinations of sigmoid functions. In discrete time, a simple piecewise affinc 
model can be defined by the following relation 

= ax\ + (1 - a) K ij H (siM ~ !«)), i=l,N (1) 

i.e. x t+1 = F(x t ) where x 1 = {%l}i=i,N is the local variable vector at time t £ Z and Fi(x) = axi + (1 — 
a) KijH(sij(xj — Tij)) is the ith component of the mapping F defined from the phase space R N into 

itself. 

In expression the subscript i labels a gene (N denotes the number of genes involved in the network). The 
graph supporting the network (the arrows between genes) is implicitly given by the sets C {].,••• , N}. 
For each i the set consists of the set of genes which have an action on i. In particular, a self-interaction 
(loop) occurs when some set contains i. Examples of networks are given Figureni an d Figure El 

The property of the action from j £ to i is specified by a sign, the number Sij. An activation is associated 
with a positive sign, i.e. = +1, and an inhibition with a negative sign Sij = — 1. 

The degradation rate a £ [0, 1) is supposed to be identical for all genes. This assumption only serves to 
simplify calculations and the resulting expression of existence domain of given orbits (bifurcation values). 
However, the whole analysis in the paper does not depend on this assumption and extends immediately to 
the case where the degradation rate depends on i. 

The symbol H denotes the Heaviside function 

„. f if x < 
H W = { lifOO 

In order to comply with the assumption of cumulative interactions, the interaction intensities Kij are sup- 
posed to be positive. They are normalised as follows 

K ij = 1 « i=hN 

This normalisation is arbitrary and has no consequences on the dynamics (see section 15.21 for a discussion 
on normalisation). These properties imply that the (hyper)cube [0, 1] N is invariant and absorbs the orbit 
of every initial condition in M. N . In another words, every local variable x\ asymptotically (when t — > oo) 
belongs to the interval [0, 1] (see section l5~T|l . Ignoring the behaviour outside [0, 1] N , we may only consider the 
dynamics of initial conditions in this set. When normalised to [0, 1], the variable x\ should be interpreted 
as a ratio of gene product concentration produced by the regulatory process, rather than as a chemical 
concentration. Lastly, the parameters belong to the (open) interval (0, 1) and represent interaction 
thresholds (see section I5~^l for a discussion on interaction thresholds domains). 



Part A. Results 
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3 General properties of the dynamics 



In this section we present some dynamical properties of networks on arbitrary graphs, and more restrictivcly, 
on arbitrary circuits. Most properties arc preliminaries results which allow us to simplify the analysis of 
circuit dynamics to follow. The mathematical analysis of the results in this section is given in section [3] 

3.1 Symbolic dynamics of genetic regulatory networks 

Following a widespread technique in the theory of dynamical systems, the qualitative features (the structure) 
of a dynamical system can be described by using symbolic dynamics |21) . This consists in associating 
sequences of symbols with orbits. To that goal, a coding needs to be introduced which associates a label 
with each domain in phase space. In our case of piecewise affine mapping, the atoms are domains - bounded 
by discontinuity lines - where F is affine. They are naturally labelled by the (elementary) symbols 6ij = 
H(sij(xj — Tij)) G {0,1} involved in interactions. In particular, the symbols 1 depend on the number of 
genes, on the interaction graph and on the interaction signs. 

By evaluating the atom(s) the image by F of a given atom intersects, a symbolic graph is obtained which 
indicates the possible (one-step) transitions between symbols. As a consequence every point in phase space 
generates via its orbit, a symbolic sequence, its code, which corresponds to an infinite path in the symbolic 
graph. 

It may happen that a code is associated with two distinct points (absence of injectivity). For instance, the 
code associated with an initial condition in the immediate basin of attraction of a periodic point is the same 
as the periodic point code. Simple examples can be found for the self-activator, see section l4~Tl 

It may also happen that an infinite path in the symbolic graph does not correspond to any point in phase 
space. In the present framework of piecewise contracting mapping, this happens when there exists an atom 
whose image intersects several atoms (absence of Markov property). The simplest such example is the 
self-inhibitor, see section FQ1 

The injectivity of the coding map associated with F can be shown to hold in the attractor. (The attractor 
is the set of points which attracts all orbits in phase space, see section HTT1 for a definition.) Considering the 
attractor amounts to focusing on asymptotic dynamics. In applications, this is particularly relevant when 
transients are short. 

That distinct points in the attractor have distinct codes is a consequence of the following property Points 
in the attractor are completely determined by their code (and the parameters). The expression of their 
coordinates is a uniformly converging series - relation Q in section I5"T1 2 

In addition the relation Q also provides a criterion for a symbolic sequence to code for a point in the 
attractor. (If it does, the symbolic sequence is said to be admissible.) The criterion, the admissibility 
condition - relation J^J) - simply imposes that, for each t, the formal point x l computed by using relation 
(0} belongs to the atom labelled by the symbols 0f ■ and is therefore a genuine orbit point. 

In practise, the analysis of a genetic regulatory network consists in analysing the corresponding ad- 
missibility condition (based on the transition graph) in order to determine which symbolic sequences are 
admissible, possibly depending on parameters. Note that according to relation when intersected with 
any hyperplane a = constant, the admissibility domain of any given symbolic sequence reduces to a product 
of threshold parameter intervals. 

3.2 Lyapunov stability and robustness with respect to changes in parameters 

The assumption a < 1, which reflects self-degradation of genes implies that F is a piecewise contraction. 
Orbits of piecewise contractions are robust with respect to changes in initial conditions and to changes in 

We also use the term symbol for the concatenation )i=i,N,jel(i) °f elementary symbols. 
2 Points x° in the attractor turn out to have pre-images x~ ' for any i £ N (see Proposition 15. 21 . Therefore, the symbols 
associated with points in the attractor are defined for all t G Z and not only for all t £ R, see relation HI . 
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Figure 1: A feedback circuit of N genes. 

parameters. Such robustness have been identified as generic features in various genetic regulatory networks, 
see ^H) and references therein. 

The robustness properties concern orbits, in the attractor of F, not intersecting discontinuities (x* ^ 
for all t € N, i = 1, N and j e The first property is Lyapunov stability (robustness with respect to 

changes in initial conditions). For simplicity, let x = {5;}i=i.Ar be a fixed point not intersecting discontinuities 

and let 5 = min \xj — T;A > 0. Then x is asymptotically stable and its immediate basin of attraction 

i=i,jv, jei(i) 3 3 

contains the following cube 

{x e 1^ : \xi -Xi\<6, i = 1,N}. 

The second property is structural stability (robustness with respect to changes in parameters). For simplicity, 
assume once again that the fixed point x exists for the parameters (a, Kij, Tij) and let 9 be the corresponding 
code. 3 In other words, assume that x computed by using relation with 6, namely 

Xi = (1 - a) E E K *ih = E K *ih> (2) 
k=o jei(i) jei(i) 

satisfies the relation 8ij — H(sij(xj — T^)) (admissibility condition JSJl for constant codes). If x does 
not intersect discontinuities (6 = min \xj — T;A > 0), then for any threshold set \T'A such that 

i=l,N, jGI(i) 3 

|T/ 3 — Tij | < 8, we have §ij = H(sij (xj — T^)) and the fixed point x persists for the parameters (a, ify , ZJy). 
Moreover, the fixed point expression J5J) depends continuously on the parameters a and Kij. Therefore, 
the code 9 also satisfies the admissibility condition 8ij = H(sij(xj — T^)) for the parameters (af,K!j,T!j) 
sufficiently close to (a, Kij, Tij). In short terms, every fixed point not intersecting discontinuities can be 
continued for small perturbations of parameters. 

Both Lyapunov stability and structural stability extend to any periodic orbit not intersecting discontinu- 
ities, see section [7| for complete statements and further perturbation results. These properties depend only 
on the distance between orbits and discontinuities. Therefore, they may apply uniformly to all orbits in the 
attractor in which case the complete dynamics is robust under small perturbations. 

3.3 Circuits and their symmetries 

The simplest regulatory networks are feedback circuits whose graphs consist of periodic cycles of unidi- 
rectional interactions. In the present formalism a circuit of length N (TV-circuit) can be represented by a 
periodic network with TV genes where = {i — 1} for all i G Z/iVZ (FigurcQJ. When simplifying notations, 
in TV-circuits the relation JJJ becomes 

xl +1 =ax t i + (l-a)H(s i - 1 (x t i_ 1 -T i _ 1 )), ieZ/NZ. (3) 

Circuits are not only interesting in their own. The dynamics of arbitrary networks can be described, in some 
regions of parameters, as the combination of dynamics of independent circuits [28], see section[7]for detailed 
statements. 

3 Every fixed point belongs to the attractor and the corresponding code is a constant sequence. 
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Figure 2: A network with 3 nodes and the corresponding network obtained by flipping the signs of incoming 
and outgoing arrows from the node 1, excepted the self-interaction sign. 

A iV-circuit is specified by interaction signs {sijigz/ATZj by interaction thresholds {Ti} ig z/./vz an d by the 
degradation rate a. With as many parameters as they are, the number of cases to be investigated is large. 
However, symmetry transformations apply which allow us to considerably reduce this number. 4 We present 
separately flips of interaction signs and transformations of interaction thresholds (internal symmetries). The 
former reduce the number of interaction sign vectors to be considered whereas the latter reduce the threshold 
domains to be studied. Details of the related mathematical analysis are given section 

3.3.1 Flipping interaction signs 

The Heaviside function possesses the following (quasi-)symmetry: 

H(-x) = 1 - H(x) for all x ^ 0. 

This property has a consequence on the dynamics of arbitrary networks (and not only of circuits): The 
mapping obtained by flipping the sign Sij of every incoming and every outgoing arrow from a fixed node, 
with the exception of any self-interaction, has (almost) the same dynamics as the original model. Precisely, 
to every orbit of the original mapping not intersecting discontinuities corresponds a unique orbit of the 
new mapping (also not intersecting discontinuities). In particular, when the attractors do not intersect 
discontinuities, the asymptotic dynamics are topologically conjugated. We refer to Lcmma l5.4l for a complete 
statement and to Figure |21 for an example of two related networks. 

In circuits, this result implies that (with the exception of orbits intersecting discontinuities) the dynamics 
only depends on the product of signs (positive or negative circuit); a property which has been largely 

igZ/JVZ 

acknowledged in the literature, see e.g. . 

3.3.2 Internal circuit symmetries 

Lemma |5.4I in section ED does not only serve to match the dynamics of circuits with flipped signs. It can 
be also applied to deduce a parameter symmetry in a circuit with fixed signs. Let S be the symmetry in 
with respect to the point with all coordinates equal to \ (xj = | for all i S Z/iVZ). The image by S 
of a circuit orbit not intersecting discontinuities which exists for thresholds T = {7i}j g z/Arz * s an orbit not 
intersecting discontinuities which exists for thresholds S(T). 

Depending on the product of signs, other symmetries follow from essentially cyclic permutations. The 
map R defined by (Rx)i ~ Xi-\ is a cyclic permutation in l w . As a representative of positive iV-circuits 
we consider the iV-circuit with all interactions signs equal to 1 : The image by R of an orbit not intersecting 
discontinuities which exists for the thresholds T = {Ti} is an orbit not intersecting discontinuities which 
exists for the thresholds R(T). By repeating the argument, additional orbits {R k (x t )} (unless R k (x t ) = x t 
for some k = 1,7V — 1) can be obtained in this circuit. 

As a representative of negative iV-circuits we consider a iV-circuit with all signs equal to 1, excepted 
sjy = —1- Let a be the symmetry with respect to the hyperplane x\ = \: The image by o~oR of an orbit not 

4 By symmetries, we mean transformations acting on the original network, not on the subsequent dynamical graphs as in [S]. 




Figure 3: (a) Graph of the map F(x) = ax + (1 — a)H(x). The point is a ghost fixed point. It attracts all 
initial conditions x < but is not a fixed point. The point 1, obtained from by applying the symmetry S 
is a true fixed point, (b) Assuming H (0) = instead of H(Q) = 1 in the definition of F changes the ghost 
fixed point to a genuine fixed point. 

intersecting discontinuities which exists for thresholds T is an orbit not intersecting discontinuities which 
exists for thresholds (a o R)(T). As before, additional orbits can be obtained in this circuit (provided that 
the original orbit has low or no symmetry) by applying a o R repeatedly. 

3.4 Ghost orbits 

In the previous section, the condition of non-intersection of orbits with discontinuities is due to lack of 
symmetry of the Heaviside function at the origin (H(— 0) 7^ 1 — H(Q)). Applying a symmetry transformation 
to an orbit intersecting some discontinuities may result in a ghost orbit (and vice- versa). A ghost orbit is a 
sequence in phase space, which is not an orbit of F, but which would be an orbit of a suitable alteration of F 
on some discontinuities (alterations which consists in letting H(0) = instead of H(0) = 1). In particular, 
a ghost orbit must intersect some discontinuities. The simplest example is the ghost fixed point which 
occurs for T = in the self-activator (section l4.1|l . see Figure 01 (For every T > 0, the point is a stable 
fixed point.) 

In spite of not being orbits, ghost orbits may be relevant in applications when they attract open sets 
of initial conditions. This is the case of the ghost orbits in circuits analysed below. These are periodic 
sequences which exist for parameters in some boundaries of domains where periodic orbits exist with the 
same symbolic sequence. Consequently, they occur for exceptional values of parameters. For parameters 
inside the domains, the periodic orbits do not intersect discontinuities and are Lyapunov stable (see section 
13.2(1 . Ghost orbits inherit this stability with a restricted basin of attraction (semi- neighbourhood). 5 

3.5 Comparison with other models: time delays 

For a = 0, the model Q becomes equivalent to a dynamical system on a (finite) discrete phase space, a 
so-called " logical network" in the literature of regulatory process models ^1 . Indeed for a = the 
symbols 0jy~ , can be computed by using only the symbols and not the variables x*- themselves. This is 
because the quantities x* +1 themselves can be computed by only using the symbols 0* -. 

For a > 0, the model is no longer equivalent to a boolean network. Indeed, for a > one needs the 
knowledge of the variables x* - and not only of the symbols - in order to compute the next state . (It 
may happen however that the system restricted to its attractor is equivalent to a dynamical system on a 
finite state space.) 

In any case, the model can be viewed as a discrete time analogue of the following system of delay 
differential equations 

5 For the ghost fixed point in the self-activator, the restricted basin of attraction is the half-line (— 00, 0[, see in Figure 1^1 
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Indeed, the behaviours of (JJJ and of the system of delay differential equations (r > 0) are remarkably similar. 
This is confirmed by results of the dynamics on circuits in the section below. On the opposite, in general, 
there are qualitative differences between the behaviours of the model Q and of the corresponding differential 
equation without delay (r = 0). 

More generally, the differences in behaviour between systems with and without delays and the similarities 
in behaviour between delay differential equations and discrete dynamical systems have been thoroughly 
analysed in the literature, see e.g. for introductory textbooks [101 115j . In particular, delays are known to 
result in the occurrence of oscillations in systems which, without delays, have only stationary asymptotic 
solutions. We refer to |17j for numerical example (in agreement with experimental results) of oscillations 
induced by delay in a negative feedback genetic regulatory circuit. As far as the modelisation of a real system 
is concerned, representing the dynamics of a regulatory network by a model with delays is relevant when the 
response times of various processes involved in the regulation stages (transcription, translation, etc.) are not 
negligible (again see for a relevant example). 

In the model JU, the degradation rate a plays the role of a delay parameter: the smaller a, the stronger the 
delay is. The examples below indeed show that the oscillations (which are absent in the differential equation) 
are more likely to occur when a is small. On the opposite, when a tends to 1, the dynamics coincides with 
the dynamics of differential equation. 

4 Dynamics of simplest circuits 

In this section, results of dynamical analysis of the four simplest circuits are presented. For the 1-circuits 
and the positive 2-circuit, the analysis is complete both in phase space and in parameter space. Due to a 
rich and elaborated phenomenology, the results on the negative 2-circuit are however only partial. Still, they 
provide a description of the dynamics over large domains of phase space and of parameter space. 

4.1 The self-activator 

The simplest network is the self-activator, namely the circuit with one node and positive self-interaction. In 
this case, the mapping F becomes the one-dimensional map F(x) — ax + (1 — a)H(x — T) (where < T < 1) 
which has very simple dynamics. Either x < T and the subsequent orbit exponentially converges to the fixed 
point 0. Or x ^ T and the orbit exponentially converges to the fixed point 1. 

In terms of symbolic graph, this means that the only possible symbolic sequences are the sequence with 
all symbols equal to (resp. to 1). The first sequence corresponds to the fixed point 0, the second to 1. 
Obviously, both sequences are always admissible. 

The dynamics of the self-activator is the same as those of the corresponding boolean network and of the 
corresponding ordinary differential equation |2fi| . In all cases, the self-activator is a bistable system. In short 
terms, delays have no qualitative influence on the self-activator. In this one-dimensional system, the reason 
is that, independently of the presence of delays, every orbit stays forever in its original atom and no orbit 
crosses the discontinuity. 

4.2 The self-inhibitor 

For the self-inhibitor (the circuit with one node and negative self-interaction) the mapping becomes F(x) = 
ax + (1 — a)H(T — x) whose graph is given in Figure Its asymptotic dynamics turns out to be the 
same as that of a piecewise affine contracting rotation. Piccewise affine contracting rotations have been 
thoroughly investigated El HH ■ The results presented here are immediate consequences of those in j^] . 
The mathematical details are given in section r6.il 

As indicated Figure 0J the iterations of every initial condition in [0, T] cannot stay forever in this set 
and eventually enter in the complementary interval (T, 1] (in the sub-interval (T,aT + 1 — a] precisely). 
Conversely, every point x° G (T, 1] has an iteration x l G [0, T] (which indeed belongs to the sub-interval 
(aT, Tj). Therefore all orbits oscillate between the two (sub-)intervals. 
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Figure 4: Graph of the map F(x) = ax + (1 — a)H(T — x) together with the invariant absorbing interval 
(aT, aT + 1 — a] . 



1/(0.8, T) 



T 
(a) 



(b) 



Figure 5: The self-inhibitor, (a) Graph of the map T i— > i/(a,T) for a = 0.8. (b) Gray-level plot of the map 
(a,T) i — > v(a,T) (White = - Black = 1). 



In terms of symbolic dynamics, it means that the symbolic graph is the complete graph (that is to say and 
1 can both be followed by either or 1). However, only special paths correspond to admissible sequences. 
Indeed, the analysis of the admissibility condition in this case proves that the only admissible sequences are 
the codes generated by a rigid rotation on the circle (x i— »■ x + v mod 1) with a unique rotation number v. 

Back to the phase space, the corresponding orbits themselves are given, up to a change of variable, by 
such a rigid rotation. Moreover, they attract all initial conditions fTheorem l6.1J) . Therefore, the asymptotic 
dynamics is entirely characterised by the rotation number which corresponds to the mean fraction of iterations 
spent in the interval (T, aT + 1 — a] . 

For almost all values of parameters, the asymptotic orbits are genuine orbits. But in a set of parameters 
(a,T) with zero Lebesgue measure, all orbits approach a unique ghost periodic orbit. 6 

The rotation number v(a, T) depends continuously on a and on T (small changes in parameters induce 
small changes in the rotation number). Moreover for a > all maps T i— > v(a,T) are decreasing with range 
(0, 1) and have a peculiar structure called a Devil's staircase, see Figure (a). 

This structure combines continuity and the existence of intervals (plateaus) where the rotation number is a 
constant rational number. Plateaus are due to structural stability of orbits not intersecting discontinuities 
( section 13. 2|) . Indeed, when the rotation number is rational, the map F has a periodic orbit which persists 
under small (suitable) perturbations of parameters. Since the rotation number does not depend on the initial 

6 In this case, the rotation number is still well-defined - and is a rational number - but the attractor is empty. 
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Figure 6: Trajectories in the continuous time self-inhibitor without delay (a) and with delay r = 1 (b). 



condition, it remains constant while this periodic orbit persists and we have a plateau. 7 

Additional properties of the rotation number are the symmetry v(a, 1—T) = l—v(a, T) which is a consequence 
of the symmetry map S (section l3.3.2f> and the unique plateau v(0,T) = | for a = (Figure [SJ (b)). 
Therefore, for a = and every T G (0, 1), every orbit asymptotically approaches a unique 2-periodic orbit 
and the attractor is just the same as in the corresponding boolean model [2fij . 

As argued in section ED)! the permanent oscillations in the dynamics of F can be attributed to the presence 
of delays. This is confirmed by the oscillating behaviour of solutions of the delay differential equation 

^ = -x(t) + H(T-x(t-l)) 

see Figure|S|(b). The behaviour of the corresponding system without delay is qualitatively different. Indeed 
for any T every trajectory of the ordinary differential equation 

^ = -x(t)+H(T-x(t)) 

converges to the globally attracting stationary point x = T (see Figure El (a)). 8 

In addition, the oscillations of F belong to the absorbing interval (aT, aT+ 1 — a] and this interval reduces 
to the point T in the limit a — ► 1. Under this point of view, the attractor of F converges in the limit of 
vanishing delay to the attractor of the differential equation without delay. 



4.3 The positive 2-circuit 

According to the symmetry of flipping all interaction signs (section 13.3. l|l the positive 2-circuit can be 
obtained by choosing, in relation (EJJ with N = 2, either si = S2 = 1 or s% — s 2 = — 1. In order to compare 
with results on other models in the literature 0121], we have opted for the second choice. The system (EJJ) 
then becomes the following system of cross-inhibitions 

J x\ +1 = ax\ + (1 - a)H(T 2 - x\) 
\ xl +1 = ax\ + (1 - a)H{T x - x\) 

which can be viewed as the iterations of a mapping F of the square [0, l] 2 . For such mapping, the symbolic 
coding follows from the partition of the square into 4 atoms labelled by 00, 01, 10 and by 11, see Figure 
[3(a). The corresponding symbolic graph is given Figure 13(b). According to this graph, for every initial 
condition in [0, l] 2 , one of the following assertions holds 

7 Ghost periodic orbits occur for T at the right boundary of plateaus. 
8 Defming H(0) = T instead of H(0) = 1. 
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Figure 7: The positive 2-circuit. (a) Atoms in phase space with affine dynamics and directions of motions, 
(b) Associated symbolic graph. 




Ti 

a = 0.25 a = 0.5 a = 0.75 



Figure 8: The positive 2-circuit. Gray level plots of the domains in the threshold plane (Ti,T 2 ) where the 
diagonal contains an orbit with rotation number v (White: v = - Black: v = 1). A projection of the 
domains in the plane (a, Ti) (or in the plane (a,T 2 )) is given Figure|5] (b). 

• either the orbit visits only the atoms 00 and 11, 

• or the orbit enters one of the atoms 01 or 10 and is trapped inside it forever. 

When trapped in 01 the orbit converges to (1,0) (lower right corner of the square). When in 10 the orbit 
converges to (0, 1). The fixed points (1, 0) and (0, 1) exist and are stable for any parameter a G [0, 1) and 
Tt,T 2 G (0,1). 

That an orbit visits only 00 and 11 depends on the parameters. If this happens, the orbit must oscillate 
between the two atoms and must asymptotically approach the diagonal, i.e. x\ — x 2 — ► when t — > oo. 9 
An orbit on the diagonal is characterised by X-^ - — ■ X<2 — X • According to the definition of F, the quantity x 1 
must satisfy the dynamics of the self-inhibitor simultaneously for T = T\ and for T = Ti 

x t+1 = ax 1 + (1 - a)H(Ti - x l ) = ax 1 + (1 - a)H(T 2 - x*) t G N 

The results on the self-inhibitor imply that these equalities hold for some x° G [0, 1] iff the thresholds T\ 
and T2 are such that v(a,Ti) = v(a,T2). According to relation JHJ in section RTT1 the attractor (or a ghost 
periodic orbit) of the system of mutual inhibitions intersects the diagonal iff T± and T 2 belong to the interval 
(the point if v is irrational) [T(a, u),T(a, v — 0)] for some v G (0, 1), see Figure|Hl 

When a periodic orbit on the diagonal exists, its basin of attraction can be determined. For instance, there 
exists a domain of parameters - defined by the inequality a < min{Ti, 1 — Ti, T 2 , 1 — T 2 } - for which the 

9 This follows from the fact that if (x\,x\) £ 00 U 11, then x\ +1 - x l 2 +1 = a(x\ - x\). 
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Figure 9: The negative 2-circuit. (a) Atoms in phase space with afhne dynamics and directions of motions, 
(b) Associated symbolic graph. 

image of 00 is contained in 11 and the image of 11 is contained in 00. In other words, for such parameters, 
there exists a 2-periodic orbit on the diagonal which attracts every initial condition in the atoms 00 and 11. 

As a consequence of the continuous dependence of v(a,T) with a, the squares \T(a,v),T(a,v — 0)] 2 vary 
continuously with a. For a = 0, the central square corresponding to v = 1/2 coincides with [0, l] 2 . Just as 
in the corresponding boolean network |2(i| . the positive 2-circuit possesses for a = two stable fixed points 
and a stable 2-periodic orbit. 

These permanent oscillations do not occur in the corresponding system of coupled ordinary differential 
equation which only present two stable fixed points (and possibly a hyperbolic fixed point on a scparatrix) 
As before, oscillations can be attributed to a delay effect due to discreteness of time in the model 
l[Tj). That a time delay is necessary to obtain permanent oscillations in positive circuits has already been 
acknowledged in the literature [111 I2()| . 

In the limit a -»• 1, the union of rectangles [aTi,Ti] x [aT 2 ,T 2 ] U [Ti,aTi + 1 - a] x [T 2 ,aT 2 + 1 - a] 
which contain the oscillations reduces to the point (Xi,T2). Once again in the limit of vanishing delay, the 
attractor reduces to that of the differential equation without delay. 

4.4 The negative 2-circuit 

Up to a flip of all interaction signs, the negative 2-circuit can be obtained by choosing s\ = 1 and s 2 = — 1 
in relation J2J with N — 2. As for the positive circuit, coding follows from a partition of the unit square 
into 4 atoms 00, 01, 10 and 11, see FigureEl(a). 10 The associated symbolic graph is given Figure[5] (b). As 
this figure suggests, no orbit can stay forever in an arbitrary given atom and every orbit visits sequentially 
every atom. 

Numerical simulations indicate that in most cases of initial conditions and of parameters, this recurrence is 
regular, i.e. the orbit winds regularly around the intersection (Ti, T 2 ) of interaction thresholds. Motivated by 
these numerical results, we have characterised such regular orbits and we have accomplished the mathematical 
analysis of their existence and of their parameter dependence. This analysis is reported in section f6. 21 and 
its main results are presented in the next two sections. 

4.4.1 Balanced periodic orbits 

According to the symbolic graph, the simplest regular behaviour is a periodic orbit passing the same number 
p of consecutive steps in each atom (balanced orbit). Its code is given by (00 p OP ll p 10 p )°° or formally 

10 Note however that the correspondence between labels and atoms differs from that in the positive circuit - compare Figures 
□ (a) and El (a). 
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Figure 10: The negative 2-circuit. Projections of parameter domains of existence of balanced orbits - for 
p = 1,5 — on the plane (a,Xi) (or by symmetry on the plane (a, T2)). The representation here differs from 
Figure [S] (intersections with planes a = constant in threshold space). A representation as in Figure |H1 would 
have presented nested squares centred at {\,\) and symmetric with respect to the diagonal; the number and 
the size of squares depending on a. 

11 

9 t+p = (aoR)(9 t ) for all* 6l and 6 l = 00 for all t = 1, • • • ,p 

According to expression J3J section I5"T1 the corresponding orbit has the same symmetry. Namely, for all 
i £ Z, x t+p is the image of x* under the rotation by angle £ with centre (i, i). In particular the orbit is 
4p-periodic and its components x* lie on the boundaries of a square. 

The existence domain of an arbitrary balanced orbit have been computed explicitely, see third item in section 
16.2.31 These domain have been represented on Figure ITU1 

The product structure in threshold space and the rotation symmetry a o R imply that, when non-empty, the 
existence domain in the plane (Ti,T2) is a square centred at (5, ^), symmetric with respect to the diagonal 
and which depends on a. For p = 1 the square is non-empty for every a £ [0, 1) and fills the whole square 
(0, l) 2 for a — 0. For every p > 1, the square is non-empty iff a > a p . The critical value a p > 5, increases 
with p and converges to 1 when p — ► 00. 

For a = 0, the 4-periodic orbit is the unique orbit in the attractor for all pairs (Xi,T2). The dynamics 
is equivalent to the corresponding boolean network |26| . On the opposite, for a arbitrarily close to 1 and 
2~i = T-2 = 21 an arbitrary large number of stable balanced orbits coexist and we have multi-stability, see 
Figure ITT1 Their components tend to (i, ^) when a tends to 1. Thus, when a is close to 1 and T\=Ti = |, 
at large scales, the dynamics is as for the corresponding system of coupled differential equations (see |26j 
for an analysis of such system). The attractor is (concentrated in a neighbourhood of) the point (|, 5) 
and every orbit has a spiral trajectory toward this point (region), see Figure El Multi-stability occurs at a 
smaller scale, inside this neighbourhood of (5, |). 

4.4.2 Regular orbits 

Balanced periodic orbits are special cases of orbits winding with a regular motion around the point (Ti,T2). 
An orbit is said to wind regularly around (7i,T2) (regular orbit) if its code is generated by the orbit of 

n The symbol 6* denotes the pair (0\,0\). The map a o R denned in section 13.3.21 simply becomes the rotation by ^ with 
centre , ^) in the present case. It transforms (xi, 22) G [0, l] 2 into (1 — X2, x\). As a consequence, we have S = (a o R) 2 , i.e. 
there is indeed only one (independent) internal symmetry in the negative 2-circuit. 
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Figure 11: Colour plot of basins of attractions in the square [0.425, 0.575] 2 for the negative 2-circuit for 
a = 0.99 and T\ = = \ - obtained from a numerical simulation. The picture shows that there are only 
balanced orbits. Several immediate basins - characterised by products of intervals - clearly appear. For 
instance, sea-green squares correspond to the balanced orbit with p = 8, brown squares to the orbit with 
p = 9, thistle squares to the orbit with p = 10, etc. In the present case, the balanced 4p-periodic orbits are 
known to exist for p = 1, ■ • ■ ,14 (indeed we have 014 < 0.99). 




Figure 12: Colour plot of the basins of attraction in the square [0, l] 2 (complete phase space) for the negative 
2-circuit with parameters as in Figure ITT1 (i.e. a = 0.99 and T\=T^ = \). The picture clearly shows that 
any orbit has a spiral trajectory toward some balanced orbit located in the region delimited by the square 
- displayed in Figure ITT1 
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Figure 13: A regular orbit code is generated by a rigid rotation on the unit circle ihi + v mod 1 composed 
by 4 arcs. Each arc (specified by a colour) is associated with an atom of the partition. Here v = g and the 
generated code is (10 2 00 01 2 11)°° and is 6-periodic. 

the rotation x h- > x + v mod 1 on the unit circle composed by 4 arcs; each arc being associated with an atom 
00, 01, 10 or 11, see Figure[T31 

In short terms the code of a regular orbit is characterised by the 5-uplc (A, B, C, D, v) where A (resp. B, C, 
D) is the length of the arc associated with 10 (resp. 00, 01, 11) and v is the rotation number. In particular 
the balanced 4p-pcriodic orbit is a regular orbit for which the 4 arc lengths are equal to \ and the rotation 
number is equal to j^. 

The balanced 4p-periodic orbits exist when both thresholds T\, T% arc sufficiently close to \ - see Figure ITT71 
When this is not the case, the orbits spend more iterations in some atom(s) than in others. The simplest 
case is when the orbit spends the same number of steps per winding in 3 atoms and a different number of 
steps per winding in the fourth atom (3 arc lengths are equal). 

In order to cover larger parameter domains, we consider those regular orbits with only 2 (consecutive) equal 
arc lengths, say A = B. Moreover, these lengths and the rotation number are chosen so that the orbits 
spend one step per winding in each the corresponding atoms (10 and 00). 

Motivated by changes in dynamics with parameters, instead of considering regular orbits with "isolated" 
numbers of steps per winding in the two remaining atoms, we consider the families of regular orbits - called 
(p, p)-regular orbits - for which one of these numbers is a fixed arbitrary integer and for which the other 
number is a continuous parameter. 12 That is to say the length C and the rotation number are chosen so 
that the number of iterations spent per winding in 01 is p G N. The (average) number of iterations spent 
per winding in 11 is represented by the real number p 1. 13 

By analysing the admissibility condition of the corresponding codes, we have obtained the following results 
on the existence, uniqueness and parameter dependence of (j>, p)-orbits. (The proof is given in section IB .2 .31 
second item.) We need the unique real root of the polynomial a 3 + a 2 + a — 1 - denoted by a c - which is 
positive (a c ~ 0.544). 

Theorem 4.1 (Families of regular orbits and their parameter dependence. Simple case) The 

(p, p)- regular orbit exists iff (T\,T2) belongs to a unique rectangle Ii(a,p,p) x l2(a,p, p) which exists for every 
p 1 and p 1 provided that a S (0, a c ]. 

The boundaries of the intervals Ii(a,p,p) (i — 1,2) are strictly increasing functions of p. The boundaries of 
l2(a,p,p) tend to 1 when p tends to oo. Moreover l2(a,p,p) reduces to a point iff p is irrational. 

The intervals Ii(a,p,p) and 7i(a,p, p') intersect when p and p' are sufficiently close. On the other hand, we 
have l2(a,p,p) < l2(a,p,p') whenever p < p' and the union |^J/2(a, n, a) consists of an interval excepted a 

countable nowhere dense set (where we have a ghost regular periodic orbit instead). 

12 Whcn this number is not an integer, it is to be interpreted as a mean number of steps spent per winding in a given atom. 
13 Technically speaking, we have A = B = u, C = pu and by normalisation D = 1 — (p + 2)u which implies that p := 3. = 
i — (p + 2). We refer to section [ti.2l for more details. 
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Figure 14: Colour plots of the existence domains in the threshold plane (T\,T2) of some regular orbits for 
a = 0.52 and a = 0.68. The orbits are characterised by one iteration per winding in any two consecutive 
atoms, p iteration(s) per winding in a third atom and p iteration(s) per winding in the remaining atom. 
On each picture, the domains where all orbits pass one iteration in 10 and in 00 - the (p, p)-rcgular orbit - 
concern the right upper quadrant (red and blue). Light (resp. dark) red corresponds to small (resp. large) 
p, i.e. number of iteration(s) per winding in 01. Light (resp. dark) blue corresponds to small (resp. large) p, 
i.e. (mean) number of iteration(s) per winding in 11. The domains where the orbits pass one iteration in 00 
and in 01 concern the left up quadrant, etc. See text for further details. 

In other words, when p' > p, the existence domain in the threshold plane of the (p, /o)-regular orbit (the 
rectangle h(a,p, p') x I 2 (a,p, p')) lies (strictly) above and at the right of the existence domain of the (p, p)- 
regular orbit (the rectangle I\(a,p,p) x I 2 (a,p,p)). As a consequence given (a,p,Ti,T 2 ) the number p is 
unique. Moreover it is an increasing function of T 2l with a Devil's staircase structure, and which tends to 
oo when T 2 tends to 1. 

The expression of the rectangle boundaries are explicitcly known (see section I6.2|l . By numerically 
computing these quantities and by applying the symmetry a o R, we have obtained Figure 1141 This figure 
presents the existence domains in threshold space, for a = 0.52 and for a = 0.68 of all regular orbits which 
pass one iteration per winding in any two consecutive atoms, p iteration(s) per winding in a third atom and 
p iteration per winding in the remaining atom. 

On the first picture of FigurclTllfa = 0.52 < a c ), the central square (in white) is the existence domain of the 
balanced 4-periodic orbit. The series of rectangles above it (from light to dark blue) are existence domains of 
(1, p)-rcgular orbits. In particular, the first large rectangle above the white square is the (1, 2)-regular orbit 
existence domain. The small rectangle in between corresponds to the (1, |)-rcgular orbit (the rectangles 
in between and corresponding to other (1, /?)-regular orbits with 1 < p < 2 arc too thin to appear on the 
picture). The second large square above corresponds to the (1, 3)-regular orbit, etc. 

The series of rectangles extending upward at the right of the series corresponding to (1, p)-regular orbits, are 
existence domains of (2, /9)-regular orbits. In particular, the rectangle in light red at the right of the white 
square corresponds to the (2, l)-regular orbit. Obviously, it is symmetric to the rectangle corresponding to 
the (1, 2)-regular orbit. The series extending upward at the right of the series corresponding to (2, p)-regular 
orbits corresponds to (3, /9)-regular orbits, and so on. Series are visible up to p = 7. 

The second picture of Figure PHI (a = 0.68 > a c ) shows that when a increases beyond a c , some rectangle 
persist whereas other do not. This is confirmed by the next statement which claims that the rectangle 
corresponding to the pair (p,p), with any p sufficiently large, persist if p is large and do not persist if p is 
small. 

Proposition 4.2 Provided that p is sufficiently large (i.e. for any p larger than a critical value p a which 
depends on a € (0,1)), the results of Theorem \4-l\ extend to the (p, p) -regular orbits with arbitrary p larger 
than a critical value (which depends on a and onp). 

(Theorem states that p a = 1 and that the critical value of p equals 1 whenever a < a c .) When a > a c , 
both p a and the critical value of p are larger than 1. The critical integer p a tends to oo when a tends to 1. 
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ha = 1 and nc = 1, 20 tia = 2 and nc = 2, 6 ha = 3 and nc = 3, 4 



Figure 15: Colour plots of (n_a, 1, nc, p)-rcgular orbit existence domains in the threshold plane for a = 0.842. 
See text for details. 

In particular, for a = 0.68 (second picture of Figure IT4"j) . we have p a = 2 and for any p ^ 2 the critical value 
of p is (at most) 2. That is to say the (p, p)-regular orbit exists for any p ^ 2 and any p 2. On the other 
hand we have a = 0.68 > a^i (see Theorem 14 . 31 below for the definition of ai.i) and thus (p, p)-regular orbit 
with p — 1 only exist for isolated values of p. In particular, the balanced 4-pcriodic orbit exists. 

We also have obtained results on more general regular orbits than only those spending one iteration per 
winding in two consecutive atoms. Inspired by the previous statement, one may wonder about the existence 
of regular orbits, denoted (tia, ns, nc, p)-regular orbits, with an arbitrary integer number of iterations spent 
per winding in each of 3 atoms (say ha iterations in 10, hb iterations in 00 and nc iterations in 01) and 
an arbitrary mean number of iterations spent per winding in the fourth atom (p £ R, p 1 iterations on 
average in ll). 14 (Under this notation, the previous (p, p)-regular orbits are (1, l,p, p)-regular orbits). As 
claimed in the next statement, it turns out that the number p can be arbitrary large only if the number of 
iterations in the opposite atom is equal to 1. 

Theorem 4.3 (Families of regular orbits and their parameter dependence. General case) Let 

nA ?S 1 and nc ^ 1 be arbitrary integers. 

The (n a, ns, nc, p) -regular orbit can exist - upon a suitable choice of the parameters (a,T\,T2) - for any p 
in an interval of the form (p c ,oo) only if Hb = 1- 

The (nA, 1, nc, p) -regular orbit exists iff(Ti, T2) belongs to a unique rectangle I\(a, nA, nc, p)xl2(a, nA, nc, p) 
which exists provided that a € [a„ A nc ^ n A,nc\ an d that p is larger than a critical value (say p ^ p a ,n A .n c )- 
The numbers a nA n anda nA . nc are known explicitely. 

The dependence of intervals Ii(a,riA,n,c, p) (i = 1,2) on p is just as in Theorem \4-l\ 

For the proof and for explicit expressions, see section T6. 2. 31 first item. Results of the numerical computation 
of (tia, 1, nc, p)-regular orbit (and their symmetric) existence domains for a = 0.842, for 3 values of ua 
and for the values of nc such that a £ [a„ Ai „ c , an A ,n c ] are presented on Figure [TBI On these pictures 
the (tia, 1, nc, p)-regular orbit existence domains are in the right upper quadrants. The other domains 
corresponds to orbits obtained by applying the symmetry a o R, 

On the first picture of Figure ITS! we have represented domain for nA = 1 and for nc running from 1 to 20. 
Explicit calculations show that a = 0.842 > 3i, no for nc = 1,2 and 3. No continuum of domains but only 
isolated domains exist for nc = 1 and 2. For nc = 3 we have a continuous series of domains for p inside 
a finite interval (first series extending upward in the right upper quadrant). On the other hand we have 
a = 0.842 £ [eij nc ,ai,n c ] f° r au n c 4. Thus for such nc, we have a continuous series of domains for all 
p larger than a critical value (unbounded intervals). Each series for nc + 1 stands at the right of the series 
for nc- 

14 Technically speaking, we have A = n^^, B = ngu, C = tiqv and p := = — — (tia + tib + nc)- 
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On the second picture of Figure ITK1 we have ua — 2 and nc runs from 2 to 6. The parameter a — 0.842 £ 
[a 2 nc , fl2,n c ] f° r n c = 2, 3, 4, 5 and 6. The continuous series of domains for nc = 6 is quite small. 

On the third picture of Figure IT51 we have ha = 3 and nc = 3,4. The parameter a = 0.842 £ [a 3 ^ nc , a3 inc ] 
for nc = 3 and 4. The continuous series for nc = 4 is barely visible. 

Families of regular orbits with arbitrary iterations per winding in the four atoms 

The existence and the parameter dependence of arbitrary families of (nA, nB, nc, p)-regular orbits (nA > 
l,ns > 1 and nc > 1) where p is varying in an interval (p\, p 2 ) (with p 2 necessarily finite) remain to be 
investigated. According to the analysis developped in section l6~2l the corresponding existence conditions 
are explicitely known (see Proposition I6.3JI . Moreover the p-dependence on threshold parameters is as in 
Theorems 14.11 and 14.31 (see Proposition ^. 4[1 . Thus, given nA,ns,nc, pi and p 2 , one only has to check the 
values of a for which the existence conditions hold. 

In the particular case where nA = Ub = nc = P, we checked the conditions numerically for several values 
of p. A statement similar to Theorem 14 . 31 results : There exists an interval of values of a (which depends on 
p) inside which the {p,p,p, p)- regular orbit exists for every p in an interval containing p. Given (Ti,T 2 ) the 
number p is unique and its dependence on (Ti,T2) is as in Theorem ^. 1\ 15 

Along the same lines, the existence of regular orbits other than {nA, fis, nc, p)-regular ones is unknown. 
For instance, we do not know if there can be regular orbits which pass p\ £ Q \ N iteration(s) per winding 
in some atom and p2 £ Q \ N iteration(s) per winding in another atom (regular orbits with non-integer 
repetitions in two atoms). 

Non-regular orbits 

In addition to regular orbits, the negative 2-circuit may have orbits in the attractor for which the code 
cannot be interpreted as given by a rigid rotation on a circle composed of 4 arcs (non regular orbits). For 
instance for (a,Ti,T 2 ) = (0.68,0.7,0.5) (a point not filled by the domains of existence of (p, p)-regular 
orbits), there exists a periodic orbit with code (10 00 2 01 3 ll 2 10 00 01 2 11 2 )°° which has period 14. For 
(a,T u T 2 ) = (0.1,0.900995,0.9005) and for (a,T 1 ,T 2 ) = (0.6,0.58,0.5), the (simplest non regular) periodic 
code (10 00 01 11 10 00 01 2 ll 2 ), which has period 10, is admissible. 

Numerical simulations show however that the typical situation when T\ ^ T 2 is a unique orbit in the 
attractor, although there are cases with several orbits (see e.g. results on balanced orbits, especially Figure 

ED). 



Part B. Mathematical analysis 



5 Analysis of general properties 

5.1 The attractor, the global orbits and the admissibility condition 

From a mathematical point of view, relation Q is interpreted as a discrete time dynamical system in 
generated by a piecewise affine contracting map F . The object of primary interested in a dissipativc 
dynamical system is the attractor. The attractor of (M. N ,F), say A, is the largest (forward) invariant set 16 
for which there exists a bounded neighbourhood U D A so that 

+oo 

A=f]F t (U), 
(see £P for a discussion on various definitions of attractor). 

15 We have checked the existence conditions of families of (p, p,p, p)-regular orbit for p = 2,3,4 and 5. (The existence 
conditions for p = 1 are given in Theorem l4.ll ) For p G {2, 3, 4, 5}, the interval of values of a has the form (ai (p), a,2 (p)) with 
< ai(p) < U2(p) < 1. Notice that for p = p, the (p, p, p, p)-regular orbit is the p-periodic balanced orbit. 

16 In all the paper, invariant always means forward invariant, a set S is said to be invariant if F(S) C S. 
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Our first statement, which is a consequence of positivity and normalisation of the interaction weights 
Kij, shows that every orbit asymptotically approaches the cube [0, 1]^. 

Proposition 5.1 A C [0, l] N . 

Proof: In order to prove the inclusion, it suffices to show that [0, 1] N is absorbing, i.e. that the image of 
any ball of radius 5 around [0, 1] N is included in any smaller ball after a (sufficiently large) finite number of 
iterations. 

The conditions Kij ^ and Kij = 1 imply the inequalities ax* ^ x* +1 ^ ax\ + 1 — a for any x* £ M. N . 

(In particular, the inclusion ^([0,1]^) C [0,1]^ follows.) These inequalities imply that d(x* +1 , [0, 1] N ) ^ 
ad(x t , [0, 1]^) where the distance d(-, •) is induced by the norm ||x|| = maxj \xi\ for any x £ R w . By induction, 
we obtain 

where B [0 ,i]n(5) = {x £ R N : d(x, [0, 1] N ) < 6}. □ 

An efficient method to describe the attractor is by using symbolic dynamics. To that goal we need to 
consider first the set of points whose orbit is global, namely G. The set G is the set of points x £ for 
which there exists a sequence {x*}t e z, called a global orbit, so that x° = x, x* +1 = F(x t ) for all t £ Z and 

SU PtGZ < +°°- 

In short terms, points in G are those points with infinite and bounded past history. Similarly as in 0] 
one proves that a sequence {x*}t 6 z is a global orbit iff its components write 

^ = (l-a)J> fc J2 K ^°'/ ■ ( 4 ) 

fc=o jei(i) 

and we have 9\j = if (sy(x*- — Tij)) for all i = 1,7V, j £ I(i) and t 6 Z. In other words every global orbit is 
entirely characterised by its code. 

Consequently, in order to ensure that the attractor can be described by using symbolic dynamics, it 
suffices to show that it coincides with G. This is the scope of the next statement. 

Proposition 5.2 The attractor A and the set G of global orbit components coincide. 

By definition the attractor of an arbitrary dynamical system satisfies F(A) C A. This statement shows that 
in our case, the inclusion is not strict, i.e. F (A) = A. 

Proof: Assume that x° £ G and let {x'}t g z be the corresponding global orbit. The relation an d positivity 
and normalisation of the -Ky's imply that x* £ [0, 1]^ C B[ 01 ]n(5) for all t ^ where 5 > is arbitrary, i.e. 

x°e f) j*(i? [0il]W (*)), 

t=0 

and then G C f|S pt ( B [o,i]^(S)) for any S > 0. 

In the proof of Proposition 15.11 we have shown that F(£?[ 01 p (S)) C -B[ .ip (S). Consequently, the set 
PlS) ^ T *(-^[o,i] JV (^)) i s invariant and is contained in the bounded ball -B[o,ip (5). By definition of A, we 
conclude that Htt^o ^*(^[o,i] N (^)) C A for every 5 > and in particular that G C A. 

In addition, let [7 be the bounded neighbourhood involved in the definition of A, i.e. A = Htt^o F 1 ^). Since 
U is bounded, there exists 5' > such that {/ C -B[o.ip(£')- Thus A C Ht^o ^* (-^[o,ip (^')) an d by using 
the previous inclusion, we conclude that A = f)fj% ^(-Bp.ip (<*')) for some s> > °- 

In order to prove that A C G, we first show that every point x £ A has a pre-image in A. The previous 
relation shows that, if x £ A, then x £ i 7 '* +1 (i3[o ! i]N(5')) for any t ^ 1. That is to say, for every t ^ 1, there 
exists 2/ t £ £?[ ip(J') such that F{F t {y t )) = x. Given t ^ 1, let z t = F t {y t ). For every f > 1, we have 
z t £F- 1 (x)n>*( J B [0!ir (5')). 
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When a > 0, every point in M. N has a finite and bounded number of pre-images by F (at most one pre-image 
for each realisation of symbols {9ij}). When a = 0, the sets F t (B^ 01 ^N (6')) (t ^ 1) themselves are finite 
with bounded cardinality. 

Therefore in both cases, there exists a pre-image z G F~ 1 (x) such that z tk = z for every k ^ where 

is a strictly increasing sequence. In other words, we have z 6 D/L^o F tk (■B[o,i] JV 0^0) an d by invariance of the 

ball _B[ ,i] N (^')> we conclude that the pre-image z of x belongs to A. 

Now by induction, for every i£A, one constructs a sequence {x'jt^o such that x° ~ x, x t+1 = F(x t ) and 
x* € A for all t ^ — 1. In other words, a; has an infinite and bounded past history. So x <G G and then A C G. 
□ 



Proposition 15.21 implies that an orbit belongs to the attractor iff its components {x*} are given by (@J 
and 9\j = H(s i j(x t j — TV,)) for all i, j and t. Equivalcntly, a symbolic sequence {0jj} te z codes an orbit in 
A iff the numbers {a;.-} computed with (@J satisfy = H(sij(xj — Tij)) for all i, j and i. Determining the 
attractor thus amounts to determining the set of symbolic sequences which satisfy this condition. It is not 
difficult to see that this admissibility condition is equivalent to the following one. 

Admissibility condition for a symbolic sequence {9 t i; j}t£z- The numbers computed by using satisfy 

sup x\ < Tij ^ inf x\ if su = +1 
fc » ■ e* =o tez : 0^=1 J 

sup J < T« < inf if Sij = -1 i = ^N,jGl{i) (5) 



tez : ej,=i 



where < means < if the corresponding bound is attained (the suprcmum is a maximum, the infimum a 
minimum) and means ^ otherwise. 



5.2 Normalisation of parameters 

In section [5] the interaction weights have been normalised Kij = 1 for every i. This simplifying 

jei(i) 

assumption can be relaxed without modifying the dynamics (up to dilations). 

Indeed, assume that the sequence {x*} is an orbit of F (not necessarily in the attractor) for some given 
parameters {K^} and {T^} (where the weights need not be normalised). Then for any vector {cti}i = i t N 
with on ^ 0, the sequence {aiX*} is an orbit of F with parameters {ctiKij} and {ctjTij}. By choosing 
a i = ^j- — 1 K . (which is always possible because 2jej(») K*j > for every i), this new orbit becomes an 

orbit of a mapping F with the weights satisfying Kij = 1 for every i. 



Another assumption in section|2is that the interaction thresholds Ty all belong to (0, 1). This assumption 
is justified by the following result on the dynamics of circuits. 

Lemma 5.3 Assume that the network is a N -circuit. If either Tj < or Ti > 1 for some i £ Z/7VZ, then 
the attractor of F consists of a unique fixed point. 

Proof: Let i be such that Ti [0,1]. We have shown in the proof of Proposition 15.11 that the distance 
between x l and [0, l] N goes to when t increases. A simple reasoning proves that the symbol associated 
with x\ does not depend on t provided that t is sufficiently large (i.e. there exists 9i € {0, 1} such that, for 
any initial condition a; , there exists t' £ N such that H(si(x\ — Ti)) = 9i for all t ^ £'). According to the 
relation J2J), we have for all t ^ t' 

xlH =a4u + (l-a)9i, 

which implies that, for any initial condition, x\ +1 converges monotonically to 9i. 

Therefore the symbol H(si + i(x t i+1 — Tj+i)) remains constant, say equals when t is sufficiently large. 
As before this implies that x\ +2 converges monotonically to 9i + \ independently on the initial condition. By 
repeating the argument, one easily proves that every initial condition converges to a unique fixed point. □ 
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5.3 Symmetries 

As announced in Section 13.31 equivalences of dynamics between distinct networks, and between distinct 
parameter values within a given network, occur which allow to reduce the number of situations to be 
analysed. These symmetries follow essentially from the next statement. 



Lemma 5.4 Let F be a mapping given by an interaction graph with N genes, and by parameters {sij}, {Kij} 
and {Tij}, and let k be fixed. The sequence {cc*} te N is an orbit of F with x* ^ Ty- for every i,j,t iff the 
sequence {x*} tS N defined by 



x\ if i 7^ k 
L — x k if i = k 

is an orbit of F with x*- ^ Ty for every i,j,t where the parameters of F are defined by 

s ij if 3 £ an d 3 ' k j • / l an ^ - _ f — s kj tf 3 G I(k) and j ^ k 
-Sik if k G and j = k % an ^ 1 Sfefe if k £ I (k) and j = k 



by = Kij for every i,j and by 

Ti 



T i3 if j G and j ^ k 
1 — Tik if k G and j = k 



As indicated in Section l3.3l an example of a network corresponding to F is given in Figure [2J 

Proof: We check that {x'jtgN satisfies the induction induced by F. The relation H(s(x — T)) = H{— s((l — 
x) — (1 — T))) implies that for any i ^ k such that k G /(i), we have H(sik(x t k — Tik)) — H(sik(x k — Tik)). 
Hence for any i ^ k, we have 

x? 1 = ax\ + (1 -a) KiiHfrM-Tii))- 

Moreover the weight normalisation and the relation 1 — H(s(x — T)) = H(—s(x — T)) which holds for all 
x T imply 

x[ +1 = axl + (l-a) KkjH{s kj {x) - T kj )). 
jei(k) 

The Lemma is proved. □ 
When the network is a circuit, and by simplifying notations as in section 13.31 the parameters of F write 

Si if % ^ k — 1, k y = / T i it i^k 



-Si if i = k — 1, k I 1 — Tfc if i = k 

By applying repeatedly Lemma 15.41 so as to maximise the number of genes for which both incoming and 
outgoing arrows are activations, one shows that, with the exception of orbits on discontinuities, the dynamics 
of every circuit with N genes is equivalent to 

• either the dynamics of the circuit where all interactions signs are positive (positive circuit), 

• or the dynamics of the circuit where all, excepted Sn, interactions signs are positive (negative circuit). 

In order to exhibit the internal symmetry S with respect to the centre of [0,l] w , starting from the 
parameter vectors (si, s%, ■ ■ ■ , sn) and (Ti, T^, • • • , TV), one applies Lemma 15.41 with k = 1. This results in 
a circuit with parameter vectors (— Si,S2, • • • , — sn) and (1 — Ti,T2, • • • , TV). Then applying Lemma 15.41 
with k — 2 produces a circuit with parameter (si, — S2, • • • , —sn) and (1 — Ti, 1 — T%, • • • , TV). By repeating 
the process until k = N results in a circuit with parameters (si, S2, • • • , sn) and (1 — Ti, 1 — T2, • • • , 1 — Tn) 
and the desired symmetry follows (see section l3*.3.2f) . 

Finally, the symmetry permutation R is obvious in a circuit with all signs being positive. In a circuit 
where all signs but sn are positive, one has to combine the permutation with Lemma [5.41 with i = 1 in 
order to preserve the signs. The result is the a o R symmetry where 



(aoRXx) = j 1 



Xi_i if i ^ 1 
— xn if i = 1 
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6 Negative circuit analysis 



This section presents the dynamical analysis of negative feedback circuits with 1 and 2 nodes respectively. 
As said in section 14.21 the dynamics of the self- inhibitor essentially follows from results on piecewise affinc 
contracting rotations. On the other hand the analysis and results of the negative 2-circuit are fully original. 

The analysis of the positive circuit with 2 nodes relies on results on the self-activator and is left to the 
reader. (The analysis of the self-activator is trivial.) 

6.1 The piecewise affine contracting rotation 

The simplest circuit which requires a proper analysis is the circuit with one node and negative self-interaction. 
The corresponding map writes F(x) = ax + (1 — a)H(T — x) whose asymptotic dynamics takes place in the 
invariant absorbing interval (aT, aT + 1 — a] (see Figure . In order to investigate this dynamics in a 
parameter independent interval, we consider the map F defined on (0, 1] bylofol' 1 where the map 
X(x) := maps (aT, aT + 1 - a] onto (0, 1]. 

The map F can be viewed as a piecewise affine contracting rotation on the circle. Adapting the analysis 
developed in [3] for such rotations, its dynamics can be entirely described. As announced in section T4. 21 the 
basic result is that all orbits are asymptotically given by a rotation. 

Theorem 6.1 Independently of the parameters (a,T), for any initial condition x G (0, 1], we have 

lim (f\x) - 4>{yt + a)) = 

t— >+oo V / 

where a G R depends on x and where the rotation number v G [0, 1) and the function <fi : R — > [0, 1] only 
depend on the parameters but not on x. The function (f> is 1-periodic and its restriction to (0, 1] is left 
continuous and increasing. 

In particular, when the rotation number is rational, every orbit is asymptotically periodic. When the rotation 
number is irrational, every orbit is asymptotically quasi-periodic. 

The dependence of the rotation number (a, T) i— > v(a, T) on parameters is known explicitely (Hj ■ It can 
be expressed as follows 

v{a,T) = v iff T(a,v) < T < T(a,u-0), (6) 

where 17 

T(a, v) = 1 - 1 ~ Qj V a k \u(k + 1)1 . 

fe=0 

Indeed, independently of a, this expression shows that the map v i— > T(a, v) is right continuous and strictly 
decreasing. Hence, given (a,T) the previous inequalities actually define a unique number u(a, T). 

Moreover and independently of a, the map T h- > v(a, T) is decreasing and continuous. The properties 
T(a, 1) = and T(a, — 0) = 1 imply that its range is [0, 1). In addition since T(a, v) < T(a, v — 0) when v 
is rational, it has a Devil's staircase structure (see Figure EJ). 

Although the map <fi also depends on parameters, it remains unchanged when T moves in the interval 
[T(a, v), T(a, v — 0)] (given a and v fixed). 

The attracting sequences \jf>{y(a, T)t + a)}t£ti (where a is arbitrary) are orbits of F in most cases. However, 
if the interaction threshold T belongs to the right boundary of an interval with rational velocity, i.e. if 
T = T(a, v — 0) for some fgQ, all sequences {cj>(v(a, T)t + a)} tg N are translated of the same periodic ghost 
orbit. In that case, the attractor A is empty 0. The set {T(a, v — 0) : v G Q} is countable and nowhere 
dense as can be deduced from the properties of the function T(a, v) stated above. 

In all other cases, i.e. for the interaction threshold not at the right boundary of an interval with rational 
velocity, the asymptotic dynamics is semi-conjugated to the rotation R u ( a ,T){x) = x + v(a, T) mod 1 on the 
unit circle. Strictly speaking we have F o <f> = <f> o R v ^ a ,T)- 

17 [ieJ is the largest integer not larger than x. 
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6.2 The negative 2-circuit 



In the negative 2-circuit, the rotation symmetry a o R applies to orbits not intersecting discontinuities (see 
section f3. 3. 2JI . In order to extend this symmetry to all orbits (and subsequently to simplify the analysis of 
their admissibility), we consider the map F defined by 

F(xi,x 2 ) = (axi + (1 - a)9 2 ,ax 2 + (1 - a)6>i) 

where 



(1,0) 


if 


X\ 


>Ti 


and 


X2 


>T 2 


(0,0) 


if 


Xl 


^ T\ 


and 


■1*2 


>T 2 


(0,1) 


if 


X\ 


<Ti 


and 


X2 


^T 2 


(1,1) 


if 


X\ 


^ T\ 


and 


X2 


<T 2 



This map and the map induced by (JHJ) are equal everywhere except on discontinuities. They have the same 
symbolic graph (Figure EJ) and the same formal expression of orbits in the attractor (relation J3J| for the 
negative 2-circuit), namely 

oo 

x\ = (1 - a) a k e\-_\~ l for aU t G Z, i= 1, 2. (7) 

The admissibility condition associated with the present map slightly differs from (j5j) . It is entirely compatible 
with the rotation symmetry; the admissibility domain in threshold space of the image of a symbolic sequence 
by a o R is the image under the same action of the original admissibility domain. 



6.2.1 Regular codes and families of regular codes 

As reported in section T4. 41 the recurrence induced by the symbolic graph suggests to consider the symbolic 
sequences generated by the rotation x t— > x + v mod 1 [y > 0) on the unit circle composed of 4 arcs with 
length respectively A,B,C and D and corresponding to the 4 atoms of the partition (see Figure EJ- Up 
to a choice of the origin on the circle, these symbolic sequences (called regular symbolic sequences or 
regular codes) are given by 



(1,0) 


if 





^vt- 


NJ < 


A 


(0,0) 


if 


A 


^ist- 


L^J < 


A + B 


(0,1) 


if 


A + B 


^ist- 


L^J < 


A + B + C 


(1,1) 


if A - 


VB + C 


s^vt- 


L^J < 


A+ B +C + D = 1 



for all f G Z (8) 



i.e. 

6\ = 1 + \yt - (A + B + C)\ - \ vt - A\ and 6\ = l+[vt-{A + B)\ - [vt\ for all t G Z (9) 
and are denoted by the 5-uple (A, B, C, D, v). 

The parameters A, B, C, D and v must satisfy the following conditions. They are positive real numbers. 
The lengths satisfy the circle normalisation A + B + C + D = 1. Moreover, in order to generate a symbolic 
sequence compatible with the symbolic graph of Figure EO the orbit of the rotation x ^ x + v mod 1 must 
pass each arc. Hence, we must have 

v < I < 1 - 3^ for all / G {A B, C, D}, 

and 

2v < I + J ^ 1 - 2u for all I ^ J G {A, B, C, D} (10) 
which in particular imply u ^ -g. 
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By analogy with the self-inhibitor, 18 we will focus on the admissibility of the regular symbolic sequences 
for which the arc lengths depend on v in the following way. Given ua, ns, nc, nn G Z, we have 

I = njv- \n lV \ for I = A,B, C, and D (11) 

Such codes are denoted by (ua, "b, nc, nn, v). Similarly as for the self- inhibitor, the rotation number v is 
allowed to vary in some interval. The novelty is that this interval depends on the numbers ha, tib, nc, nn 
because the rotation number must be chosen so that the orbit passes every arc. 

Hence, any 4-uple (nA,ns,nc , n n) generates a family of (tia, ns, nc, nn, i^)-codes with the parameter v 
varying in an appropriate interval. Moreover the normalisation A + B + C + D = 1 holds for any code in 
the family iff we have 

tia + tib + nc + Tin = and \jiav\ + \tibv\ + \_ n c^\ + [n D v\ = — 1, (12) 

for any v in the corresponding interval. In the sequel, we always assume that these conditions hold. We 
shall regard the first condition as determining nn when tia, ns and nc have been given. 

Before entering the admissibility analysis, we notice that the families (nA,nB,nc,nn,v) contain all 
regular symbolic sequences with rational rotation number. Indeed, given a sequence (A,B,C,D,v) with v 
rational, the 4 lengths can be modified (without affecting the code) so as to satisfy / = niv — \ niv\ . In 
addition for any irrational rotation number, the sequences (tia, tib, nc. Tin, v) form a dense subset of regular 
symbolic sequences with that rotation number. Therefore the sequences (nA,nB,nc,nn,v) are generic 
regular symbolic sequences. 

The decomposition of arc lengths into / = TijV — \riiv\ however may not be unique. Indeed two distinct 
5-uples (tia, tib, tic, Tin, v) may generate the same regular symbolic sequence. However if (nA,ns,nc,nn) ^ 
(n' A ,n' B ,n c ,n D ) there exists v (arbitrarily small) such that the symbolic sequences arc (nA,nB,nc,nn,v) 
and (n' A ,n' B ,n' c ,n' D ,v) are distinct. 19 

For the sake of simplicity, we shall only consider those families of {nA,ns,nc,nn,v) for which the 
condition \njv\ = holds for 3 elements in {A, B,C, D}, say A, B and C, for all v in the corresponding 
interval. In the rest of the paper we shall assume that [ n A^J = \nsv\ = \ncv\ = for all v G 
(0) nA+rl g +nc+ i ] (the corresponding interval of v in this case). In phase space, it means that the orbit 
passes (independently of v) the same number of iterations per winding in 3 atoms, precisely nA iterations 
in 10, ns in 11 and nc in 01. 



6.2.2 Reduction of admissibility for the sequences (tia, tib, tic, Tin, v) 

In the case of the sequences (tia, tib, nc, nn,v), the admissibility condition turns out to reduce to a condition 
on a real function. To see this, given z S R, n G Z and v > 0, let 20 

oo 

<p(z, n, v) = \nv\ + (1 - a) ^ a k ( [z - (k + n + l)v\ - [z - (k + ) (13) 
which is a 1-periodic function of z. 

Given any integers n a , n@ and n 7 , let n a p = n a + np and n Q ^ 7 = n a p + n 7 . 

18 Regular symbolic sequences in the self-inhibitor: The results in section 16.11 imply that the admissible symbolic 
sequences of the one-dimensional mapping F can be viewed as being generated by the rotation x t—> x + v mod 1 on the unit 
circle composed of 2 arcs with respective length E (for the arc corresponding to the atom 1) and F = 1 — E (for the arc 
associated with the atom 0). 

Under this point of view, Theorem 16. II in that section implies that for any pair (a, T) there exists a unique v := u(a,T) such 
that the sequence generated by the rotation on the circle with arc lengths E = 1 — v and F = u, either is admissible or is the 
code of a ghost orbit. The arc lengths depend on v in the following way 

/ = nju — \niv\ for I = E and F 

where tij = —1, \n^v\ = —1, np = 1 and L n F y J = 0. 

19 The same comments hold for families of symbolic sequences in the self-inhibitor. 
20 \x~\ is the smallest integer not smaller than x. 
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Together with assumption ['t-a^J = \jibv\ = the property [A+B\ =0 implies [uabv\ = \jiav\ + [tibv] = 
and thus A + B = uabv. Similarly, one shows that B + C = ubc v i A + B + C = uabc v and, using also 
the relation (JTJJ), that B + C + D = l- n A v. 

According to the definition of F and to the expression (JHJ, the symbolic sequence (jia, tlb, tic, id, v) is 
admissible iff for all t S Z we have 

x\ > Ti and x| ^ T 2 if < z4 - [vt\ < A 

x\ ^ Ti and i| > T 2 if A < z4 - [x^tj < A + S 

ccf < T x and < T 2 if A + B ^ vt - [vt] < A + B + C 

x{ > Ti and x| < T 2 if A + B + C sC z4 - [Wj < 1 

where x\ and x\ are computed by inserting the condensed expression © of a regular code into the expression 
of global orbits. It results from relations © and that the orbit components are given by 

x\ = ipiyt — \yt\ , «ab, i 7 ) and x 2 = y(z^t — \vt\ — A, rise, v ) for all t 6 Z 

Therefore the symbolic sequence (riA,nB,n>G,ii>D,v) is admissible iff the inequalities (|14|l below holds both 
with (a, /?, 7, T) = (A, B, C, Ti) and with (a, /3, 7, T) = (B, C, D, T 2 ). The inequalities are 

tp(ist - \yt\,n a p,v) > T 
tp(vt - [vt\,n a p,v) sC T 
Lp(yt - [vt\,n a p,v) < T 
Lp(yt - [vt\,n a p,v) ^ T 

Together with the assumption [^a^J = \jibv\ = L'lcH = the next statement allows to simplify further 
the admissibility condition. 

Lemma 6.2 For \n a v\ = \npv\ = and v > 0, the inequalities {1$ are equivalent to the following ones 

<p({n a + l)v - 0,n a/3 , v) < T sC (p((n a - l)v,n a p,v) 

and 

y(n a(il v - [n a ^v\ - 0, n a p, v) < T ^ ip(n a/3l v - [n a p 7 v\ , n a p, v) 
where again < means < when the limits are attained (which happens iff v is rational) and means ^ otherwise. 



if ^5 vt — \vt\ < n a v — \n a v\ 

if n a v — \n a v\ ^ vt — \ vt\ < n a pv — \n a fjv\ 

if n a ()V — [n a /jv\ ^ vt — \ vt\ < n a/ 3^v — [n a /3-yi 

if n a ^v — \ n a ji 1 v\ ^ vt — \ vt\ < 1 



for aU t G Z (14) 



Proof: Together with v > the relations (|lt)(l and (jl 1|> and the conditions [n a v\ = [npv\ = imply that 
n a > 0, rip > and ln a pv\ = 0. When n a p > 0, the expression of ip(z, n a p, v) can be rewritten as the sum 
of two functions, namely if := tpi + (p 2 where (note that |~n a( 3z/] = 1 in the present case) 

"a/3 — 1 OO 

fi(z,n a p,v) = 1-(1 -a) ^2 a k [z-(k+l)v] and ip 2 (z, n a , 3 , v) = (l-a){a~ nc " 3 -1) ^ a fc [z-(A:+l)^J . 

k— k—n a ff 

The map 2; 1— » (z, n a/ 3, z^) is a right continuous decreasing function on [0,1) which is constant on every 
interval [jv, (j + l)v) where ^ j < n a p (note that [jv, (j + l)v) C [0, 1) for any ^ j < n a p because 
\n a pv\ = 0). It is also constant on the interval (n a pv, 1]. 

The map z t— > tp 2 (z, n a p, v) is a right continuous increasing function on [0, 1). It is a step function if v is 
rational with discontinuities at pv — \pv\ where p > n a p. When v is irrational, the discontinuities are dense 
in [0,1) and the map is strictly increasing. As a consequence, every left limit <p(z — 0,n a p,v) is attained 
when v is rational and not attained when v is irrational. 

Therefore when \ n a v\ = [npv\ = the first two lines in (|14f) are equivalent to the following inequalities 

max <p(pv — 0, n a p, v) < T < min ip(j>v^n a [},v) 
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and since n a $v < n^^v — \n a p 1 v\ (because \n a p 1 v\ — \n~ i v\ and < n^v — \n^v^), the last two lines in 
<|14|) are equivalent to the following inequalities 

ip(n af3 ~ ( v - \ n af3l v\ - 0, n a p, v) < T < f(n a ^v - \ n af)l v\ , n a p, v). 

It remains to show that 



max <p(pv — 0, n a p, v) = tp((n a + l)v — 0, n a ^, v) and min (p{pv, n a p, v) = ip((n a — n a p, v). 

We prove the second assertion; the first one follows similarly. In the case where n > and nv G (0, 1], the 
expression of the function ip becomes 



ip(z, n, v) = < 



H-oc 



,LSJ 



(a-' n -l)J2 al " i if Os^zsC 



3=1 

1 - (a- n - 1)^V^J if nv^z<l 

3=0 



(15) 



By assumptions on n a ,np and v, we have \pv\ = and \(n a p — p)v\ = 1 for every ^ p < n a . It results 
that 

ip((p + n aj3 , v) - ip{pv, n af3 , v) = (1 - a)a p ((a - ™ " 3 - l)V'(^) - l) 
for every ^ p < n a , where the function ip is defined in relation l|16|l below. 

Together with the condition ["a/J^J = 0, the constraint (|10fl implies that v - 1 2 and we have 



tp((p + l)v, n a p, v) - tp{pv, n a p, v) < (1 - a)a p ( (a - l)ip( — — ) - 1 ) < 



la/3 



since ib{ — - 



3+2 ; ~~ l-o n "/3 



— . The second assertion is proved. 



□ 



By computing the quantities of Lemma l6.2l both with (a,(3,j) = (A,B,C) and with (a,/?, 7) = (B,C,D), 
one can conclude about the existence of a product of intervals in the threshold plane (Ti, T2) inside which a 
given sequence (nA,nB,nc,rir),is) is admissible. Such existence may depend on a. For an example of such 
a computation, see the analysis of balanced orbit existence domains in the next section. 

By analogy with the self-inhibitor it is also interesting to determine those triples (jia, tib, nc) ~ an d those 
values of a - for which every sequence (ua, ns,nc, n-D, v) with v G {v\, C (0, nABa+ \ \ is admissible upon 
the choice of (Ti,T2). Indeed, in addition to admissibility, one may also be interested in describing changes 
of the rotation number with threshold parameters. 

To that goal the following function ip is useful 

00 

3=1 

This function is strictly increasing, left continuous and with discontinuities for every rational number (ip is 
continuous at every irrational number). Moreover, tp(0 + 0) = and the relation — 0J = \x~\ — 1 implies 

-f OG 

^(i/ + o) = JVfi- 1 . 

3=1 

The conditions on parameters, which ensure that any symbolic sequence in {(ha, tib, nc, no, v) : v G 
{vi,V2\] is admissible provided that {T\,T%) is suitably chosen, are given the following statement. 



Proposition 6.3 Let the integers nA,ns,nc, the degradation rate a G (0,1) and the rotation number in- 
terval (vi,v 2 ] C (0, nA s C+ i ] be given. 

There exists a set S in threshold space such that, for every v G (vi,V2], the sequence (jiA,riB,nc,nYj,v) is 
admissible, upon the choice of (Ti,T2) G S, iff the following conditions hold simultaneously 



2G 



(CI) (a- riAB - 1)(1 + a)4)(v 2 ) ^ 1 and (a~" BC - 1)(1 + a)ip(u 2 ) 1. 

(C2) < 1- a"- 4 -a" c (l - a nAB ) + a" c (l - a nAB ){a- nBC+1 - 1)^(^1 +0). 

fC3j a- nABC (l-a" BC )(l-a n ^ +1 )V'(>2) 1 - a" B . 

fQJ (a"" AB - 1)(1 - a nBO+1 )ip(v 2 ) < 1 - a-' lA+1 + a nBC (aT nAB - 1). 

(C5) < 1 - a-" B+1 + (a-" BC - l)( a -"- 4B+1 - + 0). 

Proof: It follows from the previous section that there exists (Ti,T2) such that the regular symbolic sequence 
tib, tic, n Di v) is admissible iff the following conditions hold both with (a, /3, 7) = (vl, B,C) and with 
{a,p,i) = {B,C,D) 

< Ai(n a ,np, v) := (p((n a — l)f, n a p, v) — (p((n a + — 0, n a| a, 

< A 2 (n Q ,n /3 ,n 7 ,j/) := ip(n aSil v - L"-a/3 7 ^J , "-q/3, ^) - v((n Q + \)v - 0, n aj g, 1/) (17) 
< A 3 (n Q ,n /3 ,n 7 ,i/) := ip((n a - l)v,n a p,v) - ip(n a/3l v - [n a p y v\ - 0,n a p,v) 

where now the symbol < means < if v is rational and ^ if v is irrational. (The expected fourth condition, 
namely 

< ip(n a /3^v - \n a p 1 v\ , n a p, v) ~ (p(n a!3l v - \n a p 1 v\ - 0, n a p, v) 

always holds because n a p~v — \n a p 1 v\ belongs to the interval (n a pv, 1] where (p is increasing, see proof of 
Lemma 16. 21 ) 

The proof thus consists in analysing the conditions 1)17(1 with {a, (3, 7) = (A,B,C) and with (a, /3, 7) = 
(B, C, D) to obtain the desired condition. 

Analysis of the condition < Ai(n Q , np,u). Similar calculations to those in the proof of Lemma |6 . 21 show 
that 

Ax(n Q , rip, v) = a™"" 1 (l - a - {a' 71 ^ - l)(tp(u) - a 2 ^{v + 0))) 
and by left continuity of tp 

Ai(n a , np, v - 0) = a""" 1 ^ - a) (l - (a"""* 3 - 1)(1 + o)^»(i/)) 

Since n a( g > 1, we have a"" ' — 1 > 0. Moreover \j){v) < ?/>(z/+0) for all ^ and Ai(n Q , n/3, 1/ — 0) is decreasing. 
It results that the condition ^ Ax(n a ,np,v — 0) for v = v 2 implies (the same condition with v G (vi,vq\ 
and then) < Ai(n Q , n^, v) for all 1/ G (y\, v 2 \. 

On the other hand the condition < Ai(n a ,np,v) for all 1/ G (y\,v 2 \ implies, by left continuity of ip, 
< Ai(n a , np, v — 0) for all ^ G ^2], i-c. < Ax(n a ,np, v 2 — 0). 

The conditions in (CI) follow from < Ai{n a ,np 1 u 2 — 0) by choosing (a, (3) = {A 7 B) and (a, (3) = {B,C) 
respectively. 

Analysis of the condition < A2(n Q , np,n-y,v). One shows that 

A 2 (n a ,n^n 7 ^) = l-a"«+a™-(l-a"^)(a-"^+V(^ + 0)-V^))-{ Jj"^ 1 " "^) || = ^ 

For (a, /3, 7) = (A, B,C), we have ["-7^] = |_ n c^J = 0. Thus < A2(n a , n^, n 7 , v) for all ^ G (^1,^2] iff 

< 1 - a"- 4 - a" c (l - a nAB ) + a nc (1 - a" AB )(a-" sc+1 - l)V(^) w G (^1,^2] 
from which (C2) immediately follows using that a~ nBC+1 — 1 > 0. 

For (a, f3, 7) = (B, C, D), we have = L~ n ABC v \ = — 1 since we have assumed that v ^ nABO +i • Thus 

< A 2 (nQ,, 71^, n 7 , z/) for all ^ G (y\ 1 v 2 \ iff 

0^ l-a ,lB +a-" ASC (l - a nBC )(a nAB+1 - l)ip(v) ~iv G (^1,1*2] 

which is equivalent to the condition (C3). 
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Analysis of the condition < A3(n a , rip, n 7 , v). Again similar calculations to those in the proof of Lemma 
16. 21 show that 

A 3 (n a ,rip, n^v) = a""" 1 (l - a-"° +1 + (a" n ^ - l){a n ^ +1 ip(v + 0) - 
/ a™°^- 1 (a~"^ - 1) if L«7^J=0 
+ \ if L«7^J = - 1 

For (a, /?, 7) = (A, B, C), we have [^^J = 0. Therefore < A3(n Q , np, n 7 , ^) for all v G (1^1, ^2] iff 

0< l-a- nA+1 +a" BC (a- nAB - 1) - (aT nAB - 1)(1 - a" BO+1 )^(i/) £ (1^, i/ 2 ] 

which is equivalent to the condition (C4). 

Finally for (a,/3, 7) = (B,C,D), we have [n 7 ^J = — 1. In this case < As^q,, 7^3, n 7 , v) for all ^ E (vi,V2\ 
iff 

< 1 - a-" B+1 + (a~ nBC - l){a- nAB+1 - l)tp(v) Vv G {vx,v 2 \ 
from which the condition (C5) follows. □ 

Based on Proposition ^. 31 the following statement collects the results on admissibility of a regular symbolic 
sequence and states the dependence of its rotation number with parameters. 

Proposition 6.4 Let the positive integers nA,ns,nc, the degradation rate a G (0,1) and the rotation 
number interval {yi,V2\ C (0, nA ^ c+ i ] be given. Assume that the conditions in Provosition \6.'i\ hold. 

Given v G (yi, let I\(y) and Ii(y) be the two intervals given by Lemma \fi. S\ such that the regular symbolic 
sequence (ua, tib, Tic, Tip, v) is admissible iff (T\,T2) G I\(y) x hiy). 21 

The boundaries of I\(v) and of li(y) are strictly decreasing with v. In addition, we have Iziy) < 
whenever v' < v and the union lity) consists of an interval excepted a countable nowhere dense set. 

The intervals and intersect when v and v 1 are sufficiently close. 

In particular, given ua, nc, when defined, the rotation number is unique and is a decreasing function of 
T 2 . 

Proof: The left boundary of the interval h{v) (i = 1,2) is the maximum of the quantities involved in the 
left inequalities of Lemma 16. 21 Precisely, the left boundary of I\{v) (resp. hi^)) is obtained for (a, (3, 7) = 
(A, B,C) (resp. (a, 0, 7) = (B,C,D)). By using the expression (|15fl . these quantities write 

ip((n a + l)v - 0, n a p, v) = a n « (l - a{a" n ^ - l)ip(v + 0)) 

and 

ip(n a p~ ( v - \ n a p 1 v\ ~ 0,n a p,u) = 1 - o""' 1 (a""^ - l)(a _1 (l + ln y v\) + ^{v + 0)) 
Both quantities are right continuous strictly decreasing functions of v. 

Similarly, the right boundary of the interval Ii{v) is the minimum of the quantities involved in the right 
inequalities of Lemma 16.21 By using the expression 1)15(1 ■ these quantities write 

ip((n a - l)v,n a p,v) = a! 1 -- 1 (l - (o~ n ^ - 

and 

(p(n a p 7 v - \n a p 1 v\,n a p,v) = 1 - a n "^ (a~ n °» 3 - 1)(1 + [n 1 v\ +if)(y)) 
They are left continuous strictly decreasing functions of v. 

Dependence of I%{v) on v. The interval I-iiy) is obtained for (a, /?, 7) = (B,C,D). In this case, we have 
n a p 7 v — \n a p 1 v\ = 1 — uav and \n 1 v\ = — 1. 

21 When v is irrational, the interval !%{y) reduces to a point. 
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The proof of Proposition lH"3l shows that the condition (C3) implies the following condition: ^ A2 (rig, nc, no, 
for all v' G (1/1, 1/2]. Using the definition of A2 this condition is equivalent to 

ip((riB + - 0,n BC , v') ^ tp(l - n A v' ' ,n BC , v ') Vz/ G {v\, v 2 \ 

Then the right continuity of v' > <p((tib + — 0, rise, z/), the expressions of </>(l — ri^z/ — 0, iibc, z/) and 
of ip(l — n A v' ' , riBCi v') an d the fact that [n 7 z/J = —1. imply that for any v G (v\, v 2 ) we have 

¥'((ns+l)i y -0, riBC, v) = lim tp((nB+l)v'— 0,riBC,v') < lim n^f', n BC , v') = tp(l-n A v-Q, Ubc, v) 

Consequently, the left boundary of l 2 [y) is given by ip(l — n A v — 0, ngp, z/) for any z/ G (vi, v 2 ). Similarly 
one proves that the right boundary of I 2 (v) is given by ip(l — n A v, n B c, v) for any v G (fi, v 2 \. By using the 
expressions above of <p(l — n A v — 0, n B c, v) and ip(l — n A v, n B c, v) we conclude that, when v G (fi, 1^2) is 
rational, the interval ^(i 7 ) is given by 

l2(v) = (Z(y + o),m] 

where = tp(l — n A v,n B c, v) = 1 — a~ nA (a~ nBO — l)ip(v). 22 In particular, the strict monotonicity of ip 
implies that I 2 (f) < l 2 (y') whenever v' < v. Moreover, this expression of I 2 {v) implies that I2W) 

consists of an interval excepted a countable nowhere dense set. 

Dependence of I\(v) on v: The interval h{v) is obtained for (a,/?, 7) = (A, B,C). In this case, we have 
n a f3 7 v — \n OL p~ l v\ = n AB Qi> and \n 1 v\ = 0. 

By using the expression above of ip(n AB cv — 0,n AB , v) and of <*p(n AB cv, n ABl v) we obtain the following 
inequality 

lim ip(n ABC v' - 0, n AB , v') < tp{n ABC v, n AB , v) 

In particular, we have Lp{n AB cv' — 0, n AB , v') < tp{n AB cv, n AB , v) for all v' < v sufficiently close. 

Moreover in the proof of Proposition 16.31 we showed that the left inequality in condition (CI) is equivalent 
to < Ai(n A , n B , v — 0) for all v G (v\, v 2 ). By definition of Ai, the latter is equivalent to 

lim ip((n A + - 0, n AB , z/) < ip{{n A - \)v,n AB ,v) 
u' — *u 

which implies that tp((n A + l)i>' — 0, n AB , u') < ip((n A — n AB , v) for all v' < v sufficiently close. 

Similarly, one shows that condition (C3) and the condition (C4) of Proposition 16.31 imply respectively the 
following inequalities 

lim ip{{n A + \)v' - 0,n A B,v') < tp{n ABC v,n AB ,v) 
lim ip{n A Bcv' - 0,7UbX) < y{{n A - l)v, n AB , v) 

v' — *u 

Therefore the interval and Ii(v') must intersect when v and v' are sufficiently close. □ 

6.2.3 Examples 

In this section, we check the admissibility of the sequences (n A , n B , nc, tid,v) in 3 cases. The first case is a 
family with v G (0, v 2 \ [y 2 sufficiently close to 0); the second case is a family with v G (0, nA ^ c+l ] (complete 
family) and the third case concerns the admissibility analysis of a sequence with a specific rotation number, 
namely a balanced orbit. 

1) Admissibility for every v G (0,1*2] — Proof of Theorem 14.31 Theorem 14.31 eciuivalentlv states the 
admissibility of the sequences {n A , n B ,nc, nj),v) for every v G (0, v 2 \ (where v 2 > 0) depending on threshold 
parameters. The existence domain and their dependence on the rotation number are given by Proposition 
16.41 So we only have to check the conditions (C1)-(C5) of Proposition ^. 31 

22 Moreover if 1/2 is rational, the interval is given by ^2(^2) = (max{£(^2 + 0),ip((ng + 1)^2 — MBCi^mf^)]' 
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For v\ = 0, we have \j){v\ + 0) = 0(0 + 0) = 0. Together with the assumption ns ^ 1, the condition (C5) 
with v\ = imposes ns = 1. Assuming z^i = and ns = 1, the condition (C2) becomes a™ c ^ TZ^aTt- 

We now check the conditions involving v 2 . The conditions (CI) and (C3) requires 

• 1 1 a nAC +\l-a) 

0(z/ 2 ) <k mm{ ^_ {nA+1) _ - ^ > (a -(„ c+ i) _ l)(i + a )' (l _ a «c+i)(i _ a n A +2) > 

which holds provided that ^2 is sufficiently small because the right hand side is positive and -0(0 + 0) = 0. 
Finally, the condition (C4) imposes that 1 - a,- nA+1 + a" c+1 (a~ ( ™ A+1) - 1) > (i.e. that ^""^-i < a" c ) 
and is equivalent to 

1 - a- nA+1 + a nc+1 (a~ ( " A+1) - 1) 

^ (a-(«A+l)_l)(l_ a nc+2) 

which holds provided that v 2 is sufficiently small. Consequently, there exists a positive number v 2 such that 
the condition (C1)-(C5) hold with v\ = and that number v 2 iff 

o-a" A „ 1 - a" A 
n B = l and — — r < a nc sj 



a "A+i ^ 1 _ a «A+i 

Analysing these inequalities one concludes that, for any nA,nc 1, there exist $J a„ A „ ^ a nA ,n c < 1 
such that they hold iff a nA nc $J a ^ a„ 4)no . In particular, we have o nj 1 = a± nc = 0. Theorem 14.31 is 
proved. 

It may happen that the conditions (C1)-(C5) do no longer hold when v 2 becomes sufficiently large (in the 
allowed domain) depending on a. However, for tia = 1 (or nc = 1), the calculations in the next item show 
that the conditions (C1)-(C5) hold for all v 2 in the allowed domain provided that a is sufficiently small. 

2) Admissibility of an arbitrary rotation number. Proofs of Theorem 14. II and Proposition B~21 

As in the previous proof, since V\ = the condition (C5) implies that an interval of the form (0, n ^ + J 
can be admissible only if ns = 1. For the sake of simplicity, we proceed to the analysis in the case where 
ua = 1. The results in the case where nc = 1 can be obtained by applying symmetries. 

Proof of Theorem \4-l\ Theorem 14.11 states the existence and the parameter dependence of orbits with code 
(1, l,p, —p~ 2, v) for all v G (0, ^g]- Most of this statement is given by Proposition Here, we only have 
to obtain the conditions on a under which every interval (0, is possible. 

For 7%a = Ub = 1 and nc =pwe have uabc + 1 = P + 3. For v\ = and v 2 = —^3, the conditions (C1)-(C4) 
reduce to (using again 0(0 + 0) = and that 0(^3) = fz^ffs ) 

a 3 + a 2 + a - 1 ^ 1 

; - ^ C 



1 + a 

The condition a 3 + a 2 + a 1 implies that the left inequality holds for any p ^ 1 . Under the same condition, 
we have 

1 

a sC 

1 + a 

So the previous right inequality also holds for any p 1. Therefore, when a 3 + a 2 + a 1, every regular 
orbit (1, l,p, — p — 2, j/) with ^ G (0, is admissible provided that (T0 T 2 ) is suitably chosen. □ 

Proof of Proposition \4-'A The first part of Proposition 14. 21 is a special case of Theorem 14. 31 The second part 
relies on the fact that, when a 3 + a 2 + a — 1 > 0, then the condition (C3) with v 2 = does not hold for 
any p sufficiently large. Thus for such p, the regular orbits (1, l,p, —p — 2, u) cannot exists when v is close 
to □ 

p+3 

Example a = 0.68 (second picture of Figure IT4"|l . In this case, for any p ^ 2, the condition (C3) holds for 
v 2 = ^4 but docs not hold when v 2 is close to 

3) Existence domains of balanced orbits. The symbolic sequences associated with balanced orbits are 
regular with nA = ns = nc = p and v = As said in section [4.41 the symmetry a o R implies that, 
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when non-empty, the corresponding existence domain in the threshold plane is a square centred at (|, |) 
and symmetric with respect to the diagonal. As argued in the beginning of the proof of Proposition 16.31 
such a square exists iff 23 

< Aifop,-^), 0< A 2 {p,p,p,— ) and 0< A 3 (p,p,j5,-!-) 
4p 4p 4p 

Computing these quantities explicitely, we find that the first inequality and the second one hold for any 
p 1 and any a e [0, 1). Moreover, the third inequality is equivalent to 

aP- 1 1 
> - 



1 + a 2 P 2 

When p = 1 this inequality holds for any a. A simple analysis shows that, for every p > 1, there exists 
^ < a p < 1 such that this inequality holds iff a p < a < 1. 

When this inequality holds, Lemma fo . 21 indicates that the balanced orbit exists iff 

maxM^i - 0, 2p, j-), ^ - 0, 2p, < < min{<p(^-±, 2p, -L), 2p, ^)}, i = 1, 2 

Focusing on the right bound, one shows that the map a i— > <^(|, 2p, j^) is decreasing and tends to ^ when 

a — ► 1. Moreover, we have </?(§, 2p, ^) < <p(^r, 2p, j^) for every a when p = 1. When p > 1, this inequality 
holds when a is close to 1. When a is close to a p , the converse inequality holds and in this domain the map 
a i ► (p(P^-,2p, is increasing, see Figure 



7 Stability and structural stability of global orbits 

In complement to scction l3~2l the present section contains statements on robustness of orbits not intersecting 
discontinuities, namely on their Lyapunov and structural stabilities. Lyapunov stability guarantees that the 
orbit will be observed in an experiment with suitable initial conditions. Structural stability is also relevant 
because it asserts that the present phenomenology does not depend on the details of the model - in particular 
it does not depend on its piecewise affine characteristic. Every orbit not intersecting discontinuities will also 
be observed in discrete time models generated by smooth (close to piecewise affine) mappings. 

As an extension of the first claim in section l3~2l assume that {x*}t e N is an orbit of the system at 
positive distance from discontinuities, i.e. we assume that inf^jjgj^) \Xj — > 0. Proposition 2.1 in |16j 
implies that this orbit must be eventually periodic. A simple induction shows that this periodic orbit is 
asymptotically stable. Moreover, its immediate basin of attraction consist of the balls 

{x E R N : \xj -x l A< min \x\ - T„-|, Vj = 1, N] 
i : jei(i) 



Now, if inf tji Jg /(j) |x* — T,j | > 0, then inf t>i)J - e j^ |x* — T'^ \ > for any parameters sufficiently close to 
Tij. An orbit not intersecting discontinuities remains unaffected by small changes in threshold parameters. 

Moreover, the expression Q depends continuously on a and on Kij. As a consequence if, given (a, Kij, Tij), 
an orbit satisfies inf t)ij3 - e j^ |x* — T^\ > 0, then there exists an orbit in its neighbourhood, with the same 
symbolic sequence, for any parameters (a' , K^,T-j) sufficiently close to (a, ,Tij). In short terms, every 
orbit not intersecting discontinuities can be continued for small changes in parameters. 

This continuation properties implies a kind of modularity of regulation network dynamics, as suggested in 
28 ■ Indeed any collection of independent networks, the units (which can be circuits but also more complex 
networks) can be viewed as composing a unique large network in which the weights of interactions 
between units are zero. By applying the previous continuation argument, one concludes that, in this large 

23 Notice that A 2 (p,p,p, i) = A 2 (p,p, -3p, i) and A 3 (p,p,p, i) = A 3 (p,p, -3p, {-). 
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network, any combination of unit orbits at positive distance from discontinuities is an orbit at positive 
distance from discontinuities. Therefore, it can be continued to an orbit of the large network in which all 
weights are positive - some of the weights being small (networks composed of weakly interacting units) . 

Moreover for an orbit at positive distance from discontinuities, the mapping F (F p in the case of a p- 
periodic orbit) satisfies the assumptions of the Implicit Function Theorem |29j . As a consequence the orbit 
under consideration has a (unique) continuation for any sufficiently small smooth perturbation of F. 

On the other hand, changes of a mapping in regions of phase space where a given orbit never enters (e.g. 
close to discontinuities in the present case) does not affect the orbit. By combining this argument with the 
previous one, it results that if an orbit satisfies inf tjj jej(j) \Xj — Tjj \ > then it can be continued to any 
smooth map sufficiently close to F. 

Therefore, any of the periodic orbits in the piecewise affine model can be realised in a mapping which 
involves a non-linear degradation (instead of the linear one) and whose interactions consist of combinations 
of continuous sigmoidal interactions (instead of the Heaviside function) , for instance the Hill function x i— > 
Tm +X m (for an activation) or x i— > Tr J +a , m (for an inhibition) provided that m is sufficiently large. 



8 Concluding remarks and prospectives 

In this paper, we have introduced an original model for the dynamics of regulatory networks. Just as for 
other pre-existing models (especially logical networks and ordinary differential equations), the model consists 
of a dynamical system with several variables representing gene expression levels. The dynamics is based on 
an interaction graph whose topology reflects the underlying regulatory network structure. The interactions 
themselves are combinations of sigmoid functions. 

The specificity of the present model is to have discrete time and real (continuous) variables. In particular 
this formalism allows to connect the dynamics of logical networks with the dynamics of ordinary differential 
equations. Indeed we have shown that both dynamics can be recovered by adjusting a parameter, denoted 
by a, to its extreme values. 

This parameter has been argued to quantify interaction delays which can be attributed to finite chemical 
reaction rates. In other terms, we have shown that a simple discrete time dynamical system can provide 
qualitative and quantitative insights on delay effects on the dynamics of regulatory networks. Naturally 
these systems cannot be a substitute to more elaborated models of delay (intcgro-)diffcrential equations 
whose phase space is infinite dimensional. 

Our analysis of simplest circuits has shown that, just as for logical networks and for ordinary differential 
equations, any orbit asymptotically approaches a stable fixed point or an oscillatory orbit (stable periodic 
orbit or a quasi-periodic orbit). Fixed points exist independently of parameters in positive circuits and do 
not exist in negative circuits (see |24| for related recent results in ordinary differential equations). 

Moreover delays have an effect on the dynamics of the negative self-inhibitor, of the positive 2-circuit and 
of the negative 2-circuit. These systems have permanent oscillations (depending on initial conditions for the 
positive 2-circuit). The regular oscillations can be characterised by a rotation number (together with integer 
repetitions in atoms for the negative 2-circuit) which measures the average time spent in a given region of 
phase space (atom). 

In all cases, the dependence of the rotation number on parameters has three characteristic properties 
which we recall below. Due to finite accuracy in experiments, when predicting experimental results these 
properties may become more important than the exact knowledge of the rotation number. Once again, 
they are expected to appear in genetic regulatory networks with non-negligible interaction delays (and with 
degradation rates not (or slightly) depending on genes). 

The first characteristic property is continuity. The rotation number depends continuously on parameters. It 
means that small errors on parameters have few impact on the rotation number. 

The second property is that the rotation number is constant over intervals of parameters (plateaus). Not only 
small parameter changes have few impact on the rotation number, but they are likely to keep it unchanged. 
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Figure 16: Part of the symbolic graph associated with the negative TV-circuit. The nodes represent the 
symbol vectors {0i}i=i,N- 

The last characteristic property is monotonicity with threshold parameters. Namely, when a threshold 
parameter increases the average time spent in a given region of phase space increases (or decreases when 
the rotation number is a decreasing function of the corresponding threshold parameter). Monotonicity is a 
consequence of the monotonicity of interaction sigmoid functions. 

Regular orbits in the negative TV-circuit 

To conclude this paper, we mention that the analysis of regular orbits presented in section 16.21 can be 
extended to a negative circuit with arbitrary number of genes. As a representative of negative circuit, we 
consider the TV-circuit with all signs equal to 1, excepted sn = — 1 (see section 13. 3. 2|) . 

The corresponding symbolic graph contains the subgraph displayed Figure 1161 (which for TV = 2 reduces to 
the graph of Figure [5] (b)). This subgraph (in which any transition from one symbol to the following one 
only affects one elementary symbol) is the only closed subgraph for which no path can escape [22] • (This 
subgraph is the "attractor" of the asynchronous dynamics in j22j-) According to the principle of a minimum 
number of symbol flips during iterations, the orbits corresponding to paths in this subgraph are the most 
likely to be observed in experiments. 

As in the case TV = 2, for any TV £ N one shows that no orbit can stay forever in the same atom and we have 
a recurrence. As in section T6. 21 the simplest recurrence is a regular orbit (a balanced orbit when the sojourn 
time does not depend on the atom). 

A regular orbit of the negative TV-circuit is an orbit for which the code is generated by a rotation on the 
unit circle composed of 2TV arcs, each arc being associated with a vector symbol (atom) in the subgraph. 
By assuming that the arc lengths are given as in relation l|ll|l with appropriate normalisation conditions on 
the integers nj and l/i/^J) one can proceed to the analysis of the corresponding admissibility condition to 
conclude as in Theorem 14. II and in Theorem l4. 31 about the existence of regular orbits and their changes with 
parameters. 
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